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We show by a detailed study of the mean-field approximation, the Gaussian approximation, the 
perturbation expansion, and the field-theoretic renormalization-group analysis of a (p^ theory that 
its instability fixed points with their associated instability exponents are quite probably relevant to 
the scaling and universality behavior exhibited by the first-order phase transitions in a field-driven 
scalar (j)'^ model below its critical temperature near their instability points. Finite-time scaling and 
leading corrections to scaling are considered. We also show that the instability exponents of the 
first-order phase transitions equal those of the Yang-Lee edge singularity and employ the latter to 
improve our estimates of the former. The outcomes agree well with existent numerical results. 
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I. INTRODUCTION 

Phase transitions are of great importance in almost all 
scientific fields. Traditionally, a phase transition is said 
to be nth order if the nth derivative of the chemical po- 
tentials of the two phases involved is discontinuous at the 
transition point and all lower-order derivatives are con- 
tinuous [1]. Transitions of orders higher than one often 
show, however, singularity in their transition points, and 
nowadays one usually distinguishes only first-order or dis- 
continuous phase transitions (FOPTs) from continuous 
ones that include all the others [2]. The first recorded 
continuous phase transition was the critical point of CO2 
discovered in 1869 by Andrews [3]. The renormalization- 
group (RG) theory for it and other continuous phase 
transitions in general [4] was developed, however, more 
than one century later based on scaling and universal- 
ity in the critical phenomena exhibit near such transi- 
tions [5]. Some of its results have now been predicted 
in high precisions and are consistent with high-quality 
microgravity experiments [(>] . 

FOPTs, on the other hand, have apparently a far 
longer history and appear far more frequently [7-!)]. Take 
H2O as an example, its usual phase diagram contains 
three phases: ice, water, and vapor, which are common 
to daily life. All transitions between these phases are 
FOPTs except the isolated critical point at the end of 
the line of vaporization. The first well-known theory of 
phase transitions can be traced back to van der Waals' 
equation of state [lU], which, besides the description of 
the critical point, when combined with Maxwell's equal- 
area construction [11], predicts that there is a gas-liquid 
coexistence region, within which metastable states are 
separated from unstable states by a well-defined spinodal 
curve at which the isothermal compressibility diverges. 
Gibbs' theory of equilibrium [12] further identifies two 
distinct mechanisms of equilibration for the two kinds of 
nonequilibrium states, viz., nucleation and growth [13- 
15] versus unstable growth (or spinodal decomposition 
[16] that is often referred to binary systems), respec- 



tively. This picture of the FOPTs dominates their study 
for more than one century thereafter. Nevertheless, a 
theory of nucleation is still considered to be quite good 
at present even if its predicted nucleation rate agrees with 
that measures only to within several orders of magnitude 
[14, 1-5]; and the nonlinear theory of spinodal decom- 
position [17] has yet to be substantially improved [Ki]. 
The only modification to the picture is that the crossover 
of the two modes of transition is believed to be smooth 
rather than sharp and thus neither is the spinodal curve 
sharp if exists [18, 19]. In addition, metastable and un- 
stable states are often accompanied with hysteresis, a 
non-equilibrium and nonlinear phenomenon that is dif- 
ficult to control. Accordingly, a general theory that is 
comparable to the well-developed framework of the RG 
theory for the continuous transitions has yet to be devel- 
oped for the FOPTs. 

A challenge question along this line is, instead of case 
study, whether there exist also scaling and universality 
in FOPTs that may be used as general characteristics to 
study them. In the case of equilibrium, within the frame- 
work of the RG theory, the discontinuity of extensive vari- 
ables across an FOPT has been argued to correspond to 
the existence of a discontinuity fixed point with at least 
an eigenvalue equal to the dimensionality of the system 
[20]. Trivial scaling behavior then follows [21]. In the 
dynamic perspective, there are two main occasions that 
scaling and universality emerge. One is usually referred 
to as phase-ordering kinetics in which a system evolves 
to an equilibrium thermodynamic state which consists of 
two coexisting phases from a non-equilibrium metastable 
or unstable one-phase state which originates usually from 
rapidly quenching of a one-phase, thermal equilibrium 
state [7, 22]. In such systems, one finds that in the late 
stages of growth, the structure function scales by a time- 
dependent characteristic length scale which characterizes 
the size of the growing domains [23]. Several universality 
classes have also been established in relation to Model 
A, B, etc. that are originally defined in critical dynamics 
[24]. Naturally, the RG theory has been applied to such 
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phase ordering kinetics to understand the origin of such 
dynamic scaling [22]. However, lack of a small parameter 
analogous to the case in the critical phenomena renders 
such approaches essentially a scaling analysis [22]. As 
many practical transitions are driven by an external field 
or the temperature in such a way that the systems change 
from one phase completely to the other, it has been found 
that energy dissipations and/or hysteresis exhibit scal- 
ing with respect to the sweep rate of field which serves 
as a characteristic of the irreversible processes involved 
[25-32]. An RG theory has also been adapted success- 
fully to the dynamic scaling of hysteresis in a toy vector 
model with infinite number of vector components and 
subjecting to a varying external field, arriving at scal- 
ing forms with respect to the field sweep rate that are 
determined by a zero-temperature fixed point and agree 
excellently with numerical results [33] . A subsequent at- 
tempt to a more practical Ising model using Monte Carlo 
RG method has not yet been able to determine the fixed 
point unambiguously [34] . Recently, focusing on the dy- 
namics of a generic FOPT driven by an external field in 
the 0"* model below its usual critical point, we have shown 
by a ficld-thcorctic RG method that it is governed by an 
unexpected instability fixed point of a corresponding Lp^ 
model. Accordingly, it does exhibit a distinct scaling 
and universality behavior with corresponding instability 
exponents different from the critical ones [3-S]. Here we 
shall follow a more standard approach and study in de- 
tail a mean-field theory, a Gaussian theory, a perturba- 
tion expansion about the mean-field theory, and a field- 
theoretic RG theory of the if^ theory. We also employ 
the Yang-Lee edge singularity to improve the estimates 
of the instability exponents. The good agreement of these 
estimates with existent numerical results, together with 
the mean-field numerical and analytical outcomes, shows 
that, although in the mean-field theory, the FOPT and 
the third-order transition described by the ip^ theory fall 
into opposite ground states and appear different, they 
are governed most likely by the same instability point 
and hence instability fixed points. We also discuss how 
the imaginary fixed points are reached, which is related 
to the relevancy of the fixed points. 

The rest of the paper is organized as follows. Starting 
with a usual scalar (f)^ model below its critical temper- 
ature in the presence of an external field, we derive an 
associated ip^ model that is relevant to the FOPTs in- 
volved in Sec. II. Dynamics of the model is also defined 
there. After a brief review of the dynamic field theory in 
Sec. HI to set the stage, we then present details of the 
mean-field theory in Sec. IV, followed by the Gaussian 
theory (Sec. V), which describes fluctuations around the 
mean-field theory and shows clearly a divergent correla- 
tion length and a divergent correlation time similar to 
the critical phenomena, and the one-loop perturbation 
expansions (Sec. VI), which exhibits infrared divergences 
below the upper critical dimension of 6. Detailed ex- 
position of the RG theory to deal with the divergences 
appears in Sec. VII. We then argue in Sec. VIII that the 



large-scale behavior of the scalar Lp^ theory for FOPTs 
falls in the same universality class as the Yang-Lee edge 
singularity [36] and employ its existent exponents to two 
and three loop approximations to estimate the instabil- 
ity exponents for the FOPTs. As the infrared stable 
fixed point found for the ip^ theory is purely imaginary 
numbers, we enter into discussions in Sec. IX on concern 
about whether the fixed point describes true asymptotic 
scaling behavior or just crossover. A summary is given in 
Sec. X. Two appendices are included that sketch briefly 
the formulation of a supersymmetry dynamic action (Ap- 
pendix A) and the computations of the relevant integrals 
and expansions (Appendix B). 

II. MODEL 

A. Scalar cfi^ model 

We consider a model with a usual Ginzburg-Landau 
functional 

of a scalar order parameter (j) in the presence of an exter- 
nal field H, where is a coupling constant and is positive 
for stability and r = Ci{T — Tc) is the reduced temper- 
ature with Tc being the mean-field critical temperature 
and ci a positive constant. The total free energy is 

T = - \nZ = - In j P0cxp{--H[(^]}, (2) 

where the functional integral is over all possible config- 
urations and we have absorbed the temperature factor 
into the definition of T-L. The free energy thus obtained 
is the true one in the sense that it describes equilibrium 
properties of the system concerned. Therefore, it ought 
to be a convex function of the order parameter and thus 
possesses no metastable and unstable states. Analyti- 
cal continuation of this free energy to the metastable re- 
gion gives rise to a complex free energy, whose real part 
describes the equilibrium properties of the metastable 
states and imaginary part their lifetime [37]. The free- 
energy functional or Hamiltonian (1), on the other hand, 
is itself supposed to be a result of a constrained integra- 
tion in (2) of those fields over a spatial region of size a, 
or in terms of spatial Fourier transform, 

(/)(k) = j (ix0(x) exp(— ik • x), (3) 

those fields whose wave numbers are larger than a mo- 
mentum cutoff A, which is proportional to 1/a. This 
coarse-grained procedure leads to a functional that pos- 
sesses metastable and unstable states and is appropriate 
to describe their dynamical properties [38-40] . Note that 
we have used in Eq. (3) and shall use throughout the same 
symbol for both direct and Fourier transformed spaces. 
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B. Derived model 

We shall study the FOPTs in model (1) and thus shall 
take r < throughout the paper. It is well known then 
that there is a spatially uniform spontaneous magnetiza- 
tion M, in the terminology of magnetism, below Tc even 
in the absence of H. Accordingly, it is essential to shift 
the order parameter by M . In particular, let 



(4) 



then, 



= V{rM^l2 + .gAr/4! - BM) + 



,(5) 



where V is the volume of the system, 
1 



and 



-gM^ h = H - rM 



93 gM. 



3! 



(6) 



(7) 



We shall show later on that the long- wavelength behavior 
of the system is dominated by the leading ip^ term in (5) 
Therefore, neglecting the (/J^ term, we have a derived 
ip^ model [.■!">] 



HaM = Jdr jirt^^ ^ ^ 1 



[W^]'-h^\. (8) 



Note that a (p^ model has been noticed in an RG anal- 
ysis of a mean- field spinodal fixed point [41], and has also 
been used in a nucleation theory in systems with long- 
range interactions near the mean- field spinodal point [42] . 
However, these theories are only of mean-field nature. In 
contrast, we shall take fluctuations into account. 



C. Dynamics 

Metastability is essentially kinetic in origin. In order 
to deal with it, we consider a phenomenological dynamics 
governed by the Langevin equation 



dt 



-X 



+ C, 



(9) 



i.e.. Model A in the critical dynamics [24], with a Gaus- 
sian white noise ( satisfying 

(C(x,t)) =0, 
(C(x, t)C(x', t')) - 2AJ(x - x')<5(t - t'), (10) 

or equivalently, satisfying a local functional probability 
distribution 'Dp{C) 



(10') 



where A is a kinetic coefficient. The noise is supposed 
to mimic the effects on the order parameter of those in- 
tegrated degrees of freedom, which have short relaxation 
times. Accordingly, there exists in principle also a coarse- 
grained time scale and hence a cutoff frequency. Never- 
theless, we shall neglect this constraint in the following 
as wc shall consider sufficiently long time universal be- 
havior. 



III. DYNAMICAL FIELD THEORY 

The solution (f> of the dynamic equations (9) and (10) 
depends on C and is thus stochastic and ought to be av- 
eraged over the distribution (10'). The standard prac- 
tice is to formulate the problem in terms of a dynamical 
field theory [43-48] , because then standard field-theoretic 
techniques become applicable. We shall briefly repeat in 
this section the main steps leading to the dynamic action 
in Sec. Ill A and collect the definitions of the generat- 
ing functionals for connected Green functions and vertex 
functions as well as their relationship in Sec. IIIB. The 
formulation of the theory in a supersymmetry form is left 
to Appendix A, which serves also to indicate the super- 
symmetry origin of some of the relationship. 



A. Dynamic action 

A convenient method to preform the averages is to cal- 
culate the generating functional defined as 



Z[J] 



cxp 



dxdtJ{x, t)(f>{x, t) 



(11) 



where the angle brackets denote average over (. Then 
quantities of interest such as the magnetization and 
n-point correlation functions are readily obtainable by 
derivatives with respect to the external source J, 



Ci(x) = ((/.(x)) = 



1 



sz 



Z[0] 5J{x) 

Cn{xi,. ..,Xn)= Wxi) . . . (f){Xn)) 

1 (5"Z 



J=0 



Z[0] 5J{xi)...SJ{Xn) 



(12) 
(13) 



J=0 



respectively, where we have used x^ to denote (x„,t„). 

The condition that (j> must be the solution of Eq. (9) 
can be fulfilled by inserting in Eq. (11) a Dirac delta func- 
tion, which may in turn be represented by a functional 
Fourier transform by introducing an auxiliary response 
field 4> [49]. So, 



Z[J] 



Vp{c)v<i)V^j 



exp 



dndt 



J(t) + (j} 



dt 56 , 



(14) 
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where the integration over (f> runs along the imaginary 
axis, and given formally by 



J = det 



PC 



V(f) 



det 



(15) 



is the functional Jacobian for transforming to the inte- 
gration over (j) with det denoting the determinant. Inte- 
grating over the noise leads to 



Z[J,J] = J V(PVi})Jexp + J dxdt{J(p + J4>) 



with the action C given by 
£ = / dxdt 



(16) 
(17) 



To Eq. (16) we have inserted another source J conjugate 
to (j) such that response functions which arc average of (f> 
and (j) fields can be calculated by the method similar to 
Eq. (13). In particular, the two-point response function 
Gii{x,x') = {cbix)^ix')) is 



Gii{x,x') 



1 



Z[0,0] SJ{x)SJ{x') 



.7 = 
J = 



1 Sj^jx)) 
A SH{x') 



(18) 

by noting that an external field is equivalent to the source 
J up to A. The last equality of Eq. (18) expresses the re- 
sponse to an external field and thus is the origin for call- 
ing (f) the response field. The response field in Eq. (17) 
may be integrated out by a Gaussian integral, but the re- 
sultant functional is nonlinear and inconvenient for treat- 
ment. 

At this point, there are two different methods to con- 
tinue in tackling the Jacobian (15), giving rise there- 
fore to two different forms of action. We shall proceed 
with the usual one and leave the more exotic one to Ap- 
pendix A. 

The usual method is to note that the determinant [43, 
45, 47] 



J = exp 



6(0) A / dxdt 



5^n 



= exp (-£'), (19) 



where the Hcaviside step function defined as 

e(t) = |i' 



(20) 



has to be assigned a special value 8(0) = 1/2 for con- 
sistency in the continuum limit of a symmetrized dis- 
cretization of the Langevin equation [44, 45, 47]. One 
therefore has a dynamical field theory with a total action 
Ctot = ^ + Standard field theoretical methods can 
then be utilized to compute relevant response functions. 
In particular, Feynman rules for perturbation expansions 
may be defined (see Sec. VI below). One can then prove 



that diagrams with a closed response loop given by a con- 
traction of (f) and (f> at the same spatial-temporal point 
(i.e., both (j) and come from the same term X(j>ST-l/S(f> 
in Eq. (17)) give just XG{0)S'^n/S(j)'^ = -C and thus 
just cancel the Jacobian exactly in each order of the per- 
turbation expansions because G(0) = 0(0) (sec Eq. (75) 
below). As a consequence, one may simply ignore the 
Jacobian by choosing 6(0) = [43, 44, 47], which cor- 
responds to a forward discretization, and meanwhile, ex- 
cludes closed response loops from the perturbation ex- 
pansions. Note that in this way, in the ensuing dynamical 
field theory. 



Z[J,J] = / V(f>V4>cxp 



dxdt{J(j) + J4>) 



(21) 



with the dynamic action given by Eq. (17), causality be- 
comes automatically implemented because [43, 50] 

(0(ti) . . . <^(t„)0(t'i) . . . 0(C)) = 0, if one t[, > all U. 

(22) 

In particular. 



(23) 



B. Connected Green functions and vertex functions 

From Z[J,J], Eq. (21), or, the partition function in 
statistical mechanics as it may be regarded as defining the 
statistical weight of a configure ((/), (/>), one defines another 
generating functional W[J, J], the minus free energy, by 



W[J,J] ^ In Z[jJ] 



(24) 



for connected correlation and response functions, or the 
cumulants 



{(j){xi )... (t){Xn)(i>{x'i ,...,x'„,))c 
<5"+"'W[J, J] 



SJ{xi)...6J{xn)SJix[)...Jix'^,) 



J = 
J = 



(25) 



which are generally deviations and their moments. For 
example, the connected two-point correlation function 

G^,oi^,x') = mx)-{cf,)][<i>{x')-{m 

^C2ix,x')-{^y =Cix,x') (26) 

is the second moment of the deviation — (0). Accord- 
ingly, one may construct the correlation and response 
functions from their connected counterparts, whose Feyn- 
man diagrams are fewer in number in high orders in per- 
turbation expansions. 

The connected diagrams contain one- line reducible and 
irreducible diagrams that depend on whether they are 
connected or not when an internal line of the diagrams 
is cut. The former is products of the latter, which can 



5 



be proved to be generated by the Gibbs free energy T 
through a Legendre transformation 

r[(^), (</))] = -W[J, J]+ f d^dt + J(^)) , (27) 



where (0) and (0) are averages in the presence of J and 
J given by 



SW ~ SW 



Then, from Eqs. (27) and (28), 



ST ~ ST 
J= -TT-r, J 



and the one-line irreducible vertex function Tn'n is 

Tn'ri (^1 1 • • • 1 •^n' i "^l ■>••■■) ^^n ) — 

j"'+"r[(0),(<^)] 



(28) 



(29) 



J = 
J = 



(30) 



In particular, differentiating Eq. (29) with J and J and 
using Eqs. (25) and (26) give 



ro2(a;i,a;2) = 0, (31) 
dxrii{xi,x)Gii{x,X2) = S{xi - X2), (32) 

c?a;r2o(a:^i,a;)G'ii(a;,a:;2) ^ - J dxC^oixi, x)rii{x, X2) 

= dxC{xi,x)Tii{x,X2), (33) 



where the first equation is a result of Eq. (23) and use has 
been made of Gf^ = Gn due to Eq. (23). In addition, 
the fluctuation-dissipation theorem 



\[Gii{x,x')-Gii{x',x)] = - 



dC{x,x') 
di 



(34) 



holds relating the response function to the correlation 
function [43, 45-48, 51, 52]. Moreover, dynamic response 
functions converge to their corresponding static correla- 
tion functions at long times [43, 45, 46, 52]. In fact, these 
are consequences of the supersymmetry in the theory [45] . 



IV. MEAN-FIELD THEORY 

In this section, we study the mean-field theory of a 
field-driven FOPT and show that it is controlled by an 
instability point [35], which is the spinodal point in the 
mean-field theory. 



TABLE I. Mean-field critical and instability exponents. 
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A. If model and instability point 

The mean-field theory is the lowest order of a saddle- 
point approximation to the free energy (2). Assuming 
that the saddle point lies at a uniform order parameter 
M, we have the free energy density 

F ^ J'/V = 'H/V =^rM^ + ^gM* ~ HM. (35) 
The dynamics then reduces to 



H 



(36) 



from Eqs. (9) and (35). It is well known that there is 
a mean field critical point at r = and H = with 
its mean field critical exponents collected in Table I for 
later comparison. For r < 0, on the other hand, there is 
an equilibrium FOPT at i7 = between the two phases 
with AI = —(drj g = ±Me for each given r and g. 
As there is no fluctuation, the FOPT cannot take place 
at the equilibrium transition point at iJ = because 
there is a free-energy barrier of \F{0) — F{Me)\ = /2g 
between the two phases. Rather, it can only take place 
beyond the point where the barrier vanishes. This is in 
fact the spinodal point at which both the first and the 
second derivatives of F with respect to M vanish. 

To see the essence of this transition, let us set M ~ 
Ms + m{t) in Eq. (36), which becomes 



dm , / 1 , r 2 1 s J 

— = -A I rm + -gMsUi + —gm - h 



(37) 



where Ms is a constant, and 

T = r+ ^gMl h = H- rM,, - ^gM^ (38) 

which is just Eq. (6) at Ms- One finds a pair of the 
spinodal points which lie at 



Ts =0, hs^ 0, 



(39) 



or from Eq. (38), at Ms = ±^-2r/g and Hs = 
^{2r /2,)^J —2r / g. Accordingly, Eq. (38) can be expressed 
as 



r = ci(T-T,), h = H-Hs 



(40) 



with Ts = Tc — gM^ /2ci at cither point. Comparing with 
Eq. (36), one sees that each point at its associated tran- 
sition plays a similar role to the critical point. The only 
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difference is that one has a quadratic term in Eq. (37). 
However, in the vicinity of the spinodal point, m is small. 
This term overwhelms the cubic one and thus controls the 
dynamics. Therefore, upon neglecting the cubic term in 
Eq. (37), the dynamics in the vicinity of the spinodal 
point becomes 



dm 
Itt 



-A 



(41) 



which is governed by the derived Lp^ model, Eq. (8), 
whose corresponding mean-field free energy density is 

F3 = irm^ + ^<?M,m3 - hm, (42) 

and each spinodal point, Eq. (39), now becomes the insta- 
bility point of the state m = 0, because at these points, 
the system becomes unstable as can be seen in Fig. 1(a). 

Note that if wc neglect the quadratic term in Eq. (37), 
we return to the 0^ theory. Only the critical point is 
now replaced by one of the two spinodal points. This is 
the theory of pseudo-critical phenomena in which usual 
critical behavior emerges at the spinodal point instead of 
the critical point [53]. However, in this case, the critical 
fixed point is unstable against the neglected quadratic 
term and only crossover may be observed [53]. 



B. Mean-field static instability exponents 

As the free energy of the Lp^ model is known, one can 
then derive its associated exponents similar to the 0^ 
mean-field model. We call them instability exponents 
or even spinodal exponents as they are associated with 
the instability point or spinodal point. One will see that 
they are exactly the counterparts of the critical expo- 
nents. Accordingly, we shall use identical symbols with 
the critical ones. 

At equilibrium, Eq. (41) leads to the equation of state 



1 2 
Tin + —gMsUi = h. 



As a result. 



0, 

-2T/gMs 



{T, - T) 



T > 0, 
T < 



(43) 



(44) 



for h ~ Q. So, m changes continuously to zero with j3 = 1. 
Also, as h ^ m? at r = 0, (5 = 2. Differentiating Eq. (43) 
gives rise to 



dm 
'dh 



h=0 



{T-Ts 



r > 0, 



i-r)- 



{Ts - T)-T , T < 



(45) 



for the susceptibility. So, 7 = 7' = 1, which arc identical 
to their critical counterparts. However, the amplitude 
ratio between the susceptibility above and below r = 



^3 


T= 




r>0; / / 












(a) 



FIG. 1. (Color online) (p^ mean-field free-energy density for 
(a) different rs at /i = and (b) different /is at r = 0. (c) 
Generic hysteresis of the 0^ model in an external field H. The 
arrows depict the transition that can be described by (b). The 
stars mark the spinodal/instability points. 



is now 1 instead of 2. The free energy density for the two 
solutions is 



F. = 



0, 

2tV35'M, 



r > 0, 
r < 0. 



(46) 



Thus the specific heat is 



C 



-T 



QT2 



0, 

-ATclT/gHP^ ^{Ts-Ty 



T>Q, 

T < 0, 

(47) 

which again changes continuously from one state to the 
other with a = — 1. These instability exponents are also 
collected in Table I. 

A caveat is needed here. One sees from Eq. (44) 
and Fig. 1(a) that the instability exponents just derived 
describe the continuous transition from the state with 
m = —2T/gMs to that with m = 0. As the specific heat 
is linear in r from Eq (47), this is in fact a third-order 
phase transition in the classical classification [1]. Docs 
this transition have anything to do with the FOPT we 
attempt to study? Our answer is yes. This can be seen 
from Fig. 1(b). As mentioned, there are two transitions 
associated with two spinodal points (marked by starts 
in Fig. 1(c)) which are related by inversion symmetry. 
Without loss of generality, let us consider the one with 
positive Ms and negative Hg- The FOPT is then driven 
by increasing H in the direction opposite to as shown 
in Fig. 1(c). This corresponds to a positive gMg and h 
changes from positive to negative, as Fig. 1(b) shows. 
As long as > H > Hs and hence h = H — > 0, 
the system keeps in the metastable state with M > 0. 
This corresponds to the situation in which m resides in 
the well. As h decreases, the well becomes shallower and 
closer to the origin. Exactly at = 0, the system lies at 
the transition point and becomes unstable. When h < 0, 
the positive m state loses its stability and would fall to 
negative infinity eventually. In practice, of course, the 
system restabilizes at a negative M due to the neglected 
quartic term in the free energy when m gets large (neg- 
ative). Thus, the exponents are in fact the properties 
of the transition point at h = rather than the well at 
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h > 0. This becomes apparent in the RG theory of the 
critical phenomena and also of FOPTs as will be seen 
later on. In this theory, the exponents arc properties of 
a fixed point. Only the critical point of a real system 
converges upon renormalization to the fixed point, not 
the equilibrium phase. This indicates that the RG is in- 
sensitive to the minima which arc a finite distance away 
[54], although this argument which arose from the RG 
study of the Potts model whose Landau mean-field the- 
ory predicts an FOPT instead of a continuous transition 
[55] was argued to be doubtable [54]. Nevertheless, we 
shall show in this paper that the instability point does 
describe the FOPT. However, in the FOPT we consid- 
ered, the system falls into the unstable left side instead 
of the right well. This may be the reason why the fixed 
point we shall find below is imaginary ['i^]. 



Mean-field hysteresis exponents and their 
verification 



We now show that the mean-field instability exponents 
are indeed relevant to the FOPTs considered. 

In order to probe these exponents, we shall utilize the 
method of finite-time scaling proposed in critical phe- 
nomena [5G, 57], a method in which H is swept linearly 
through the instability point. To be specific, we consider 
again the transition depicted in Fig. 1(c) and assume that 



h = -Rt, or, H = Hs- Rt, 



(48) 



where i? is a constant. This form of driving implies 
choosing i = at /i = 0. In fact, one can start sweep- 
ing the field from anywhere sufficiently far away from 
Hs- For example, H = Hq — Rt. However, /i = al- 
ways shifts the time origin to the time at H^, namely, 
h = — i?(t — tg), which recovers Eq. (48) upon setting 
ts — 0. The linear driving imposes on the system an 
effective finite time scale which is proportional to R^^ . 
When it is shorter than the correlation time, finite-time 
scaling follows in close analogy to the occurrence of finite- 
size scaling when the size of the system is smaller than 
its correlation length. We shall show below in Sec. VII G 
that at the instability point t = and h = 0, m follows 
a finite-time scaling form 



where 



m{h, R) = R^'"''" fi-hR-^^/"''"), (49) 



rH=z + (iS/u (50) 



is the RG eigenvalue associated with R, z and v are in- 
stability exponents, and / is a scaling function. Using 
the mean-field instability exponents derived in the last 
section and those that will be derived in next section, all 
collected in Table I, one finds that Eq. (49) becomes 




FIG. 2. (Color online) Mean-field transition fields Hs — Ht 
(filled symbols) and moments Mt — Ms (opened symbols) ver- 
sus R. Full lines through filled symbols and opened symbols 
are lines of slopes 2/3 and 1/3, respectively, through the data 
points at R = 10^^. The two dashed lines have slopes 3/5 
(lower) and 1/5 (upper) for comparison. Note that the sym- 
bols have small vertical variations due to the interpolations in 
extracting the points but can never fit the dashed lines. The 
parameters chosen are r = —1, A = 1, and g = 1 (circles) and 
g = 0.1 (triangles). 



This means that at m 
field 

ht = Ht — Hs 



or M ~ Mg. the transition 



C2i?2/3 ^ R- 



and at h = 0, the transition moment 

mt(0, R) = Mt - Ms = /(0)i?i/3 ~ i?" 



(52) 



(53) 



where C2 is the root of /, i.e., /(— C2) = and nn = 
P5/vrH = 2/3 and n„i = j3/vrH ~ 1/3 are hysteresis 
exponents [35]. All these results are born out by direct 
numerical solutions of Eq. (36) shown in Fig. 2. Note 
that if, instead of the theory, the original (/)^ theory 
governed the transition, the critical exponents listed in 
Table I would give uh = 3/5 and = 1/5 [58], which 
can be ruled out from Fig. 2. This indicates that the 
theory of pseudo-critical phenomena [53] cannot describe 
at least the mean-field theory. 

In fact, the above results can even be proved analyti- 
cally. Eq. (41) is a kind of the Ricatti equations that can 
be solved analytically [59]. For the driving (48), 



_2_dy 

gXy dt ' 



s = {/g\m/2t 
transform Eq. (41) into the Airy equation 



(54) 



ds2 



sy (55) 
which is solved by the Airy functions Ai and Bi as 



m{h,R) = R^/^fi-hR- 



2/3X 



(51) 



y = C3^i(s) + CiBi{s), 



(56) 
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where C3 and C4 are constants to be determined by initial 
conditions. Consequently, 




4R c^Ai'{^g\^R/2t)+CiBi'{^gX^R/2t) 
C3Ai{^/^X^R/2t) + CiBii^g\m/2t)' 



(57) 

from Eqs. (54) and (56), where the primes indicate a 
derivative with respect to the argument. Therefore, one 
recovers Eq. (52) from Eq. (54) with C2 = st/ ■\/g)^j2 
at s = St at which m = and Eq. (53) from Eq. (57) 
at h = Rt = with /(O) = {/4j^[c3Ai'{0) + 
C4Bi'{0)]/[c3Ai{0) + c4Bi{0)]. 

These mean-field results therefore shows clearly that 
the dynamic scaling of hysteresis in the FOPT of the 
(j)'^ model is determined by the instability point of the 
derived (p^ model. 



V. GAUSSIAN THEORY 



In this section, we consider the Gaussian theory for the 
ip^ model. This theory, which can be solved analytically, 
also describes fluctuations in the state m = 0. As a 
result, the other instability exponents v, z, and rj can 
all be derived. More importantly, we shall see clearly 
that at the instability point, the correlation length and 
the correlation time both diverge similar to their critical 
counterparts. 

To this end, we use the dynamic action, Eq. (17), 
which, for "Ha in the absence of the cubic term and h, 
is 



Co 



dxdt < Lp 



dip 
'dt 



+ X(t(p- V>) 



Xip 



dkdUJ { [iuj + A(r + k^)]ip{k,uj)(p{-k, -uj) - X\ip(k,uj)\^} , (58) 



where we have written (p in place of (/) for uniformity and the inverse Fourier transform is defined as 

ip{x,t) — J (ikciaJ(p(k, w) exp(ik • X — ia;t) 



with 



dkdlJ = 



1 



(27r 



\d+l 



dkdw, 



(59) 



(60) 



d being the space dimensionality. Note that, as we used a finite volume, we should in principle replace [1 / (Ztt)**] J dk by 
a discrete sum (l/V) for consistency. However, only in the thermodynamic limit in which y — > 00 and hence the 
sum becomes the continuous integral can singularity emerges. So, we still use the integral representation. Integrating 
out both (p and ip, one finds 



W[J, J] 



dkduj 



G(k, -w) J(-k, -w) J(k, Lu) + G(k, w) J(k, w) J(-k, -uj) + C(k, w) J(-k, -w) J(k, lu) (61) 



from Eq. (24), where we have neglected some irrelevant constants, 

1 



Gik,u) 



+ A(T + k2)' 

is related to the two-point response function defined in Eq. (25) by 

Gii(ki,wi;k2,W2) = J dxidx2Gii{xi,X2)e'^U('-^'^-^y^^) = (2^)'^+i(5(ki + k2)(5(wi + W2)G(ki, wi), 
the last equality being resulted from translational invariance, and 

2A 



G(k,a;) = 2AG(k, w)G(k, -w) 



;2 + A2(T + k2)2' 



(62) 



(63) 



(64) 



is related to the two-point correlation function defined in Eq. (26) through a relation similar to Eq. (63). Equa- 
tions. (62) and (64) satisfy of course the fiuctuation-dissipation theorem, Eq. (34), 



A[G(k, Lj) - G(k, -a;)] = iLuC(k, co) 
in its Fourier transformed form. From Eq. (61), one also has 

G;52(k,^)=0, 



(65) 
(66) 
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which is consistent with Eq. (23). As a result, Eqs. (32) and (33) become 

rn(k,a;) =G(-k,-a;)-i =za; + A(T + k2), r2o(k,a;) = ^ _2A. (67) 

In fact, since 

(V3(k,cj)) = G(k,cj)J(k,w) + C(k,w)J(k,w), ((^(k, cj)) = G(k, -cj) J(kw) (68) 
from Eq. (28), solving out J and J and using Eq. (27), one finds 

r= dkdij[G{k,Ljr\^{^k, -c^))(y.(k,c^))+G(k,-c^)-i(^(k,w))(^(-k,-c^)) -2A(^(-k,-c^))(^(k,w))] , 

(69) 

I 



which recovers Eq. (67) correctly by Eq. (30). Combin- 
ing Eqs. (65) and (67), one can write the fluctuation- 
dissipation theorem in the form 

iwr2o(k,c^) ^ A[rn(k, -Lu) - rn(k,c^)], (70) 

or, 

2A 

r2o(k,w)= Imrn(k,t^), (7f) 

which is of course satisfied by the Gaussian results. 

Equation (62) implies a relaxation or correlation time 
teq(k) satisfying 

Ateq(k)-(T + k2)-i=r/i(kO, (72) 

where 

^ = _ ^-'^ (73) 

is a correlation length and 

/i(x) = -l^ (74) 
f -t- a;^ 

is a scaling function with the instability exponents v = 
1/2 and z = 2. Indeed, the temporal Fourier transform 
of Eq. (62) gives 

G(k, t) = e(i)e-^("+''')* = e(t)e-*/*°<'('') (75) 

from the residue theorem by noting that G(k, oj) has a 
pole in the negative imaginary uj axis. So, the correlation 
dies out after the time teq(k). Consequently, G(k, w) can 
be cast in a scaling form [24] 

G(k,c.)=e'-V2(ke,(w/A)r)/A (76) 



with r/ = and the scaling function 

f2{x,y) = - ^— ^, (77) 

where the extra A arises from the time unit. So can 
G(k,a;) and rii(k,a;). Note that icq(k) diverges for 
|k| and t — similar to ^. 

As pointed out above, there exist relations between dy- 
namic and static functions. One sees from Eqs. (62) and 
(32) that at w = or the long-time limit, the dynamic 
functions equal the static ones up to the kinetic coeffi- 
cient, which is true beyond the Gaussian theory. Thus 

AG(k,0)=G^'(k) = ^— , (78) 

where G'^'(k) is the static correlation function indicated 
by the superscript. Its inverse Fourier transforms at d = 
3 and at t = are proportional to exp(— y^rlkj) and 
respectively, and hence confirming again Eq. (73) 
and 77 = 0. Such kinds of relation can be used as a check 
to the resultant vertex functions in dynamics as will be 
seen in the following sections. 

Note that the Gaussian theory is in fact identical for 
both 0^ and cp'^ theory. This is the reason why the expo- 
nents derived in this section are also identical. 



VI. PERTURBATION EXPANSIONS 

Now we consider the whole model Eq. (8). After 
Fourier transform, its associated dynamic action becomes 



£ = £0+2^3 J dkdk'dUdu}''f{k,uj)ip{k' ,uj')ip{—k~' k' ,—uj — uj') — X j (ikdaJ(^(— k, — w)/i(k, w), (79) 

where Cq is given by Eq. (58) and g'^ = Ags. The standard method to deal with Eq. (79) is to set up perturbation 
expansions using Feynman diagrams. One observes that Xh plays the role of J when comparing Eq. (79) with Eq. (16). 
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Accordingly, we shall not take Xh as an element of the expansions but rather regard it as a specific value of J. On 
the basis of the Gaussian theory, let a line represent ip, a dash line (p, and a filled circle — 2A, then 



G(k,-c.) 



-2A= — - 



Accordingly, to one-loop order, 

rio(k,w) = 



rn(k,w) = 



= iuj + A(t + k^) - g'i j dk'rfw'C(k', w')G(k - k', -uj + J) 




= iuj + \{T + \<?)-g'i j dk' 



1 



r2o(k, = 



T + k'^ioj + X [2t + k'2 + (k - k')2] ' 
-(^- ==-2A-igfy" dFda7C(k',w')C(k- k',tj - 



-2A - A^g^ / dk' ^ 



2t + k'2 + (k - k') 



UJ 
l\2 



T + k'2 T + (k - k')2 ^2 + [2^ + k'2 + (k - k')2]2 



(80) 

(81) 



(82) 



(83) 



ri2(ki,wi;k2,W2) = 



.93 + 5^ / dkduj [2C(k, w)G(ki - k, w - wi)G(k2 + k, w + wa) + C(k, w)G(k3 - k, w - w3)G(k - ks, ^2 - w)] 
,^ 1 



g'a + g'i J dk 



r + k2 [ iwi + A [2t + k2 + (ki - k)2] -iuj2 + A [2t + k2 + (k2 + k)2] 



ILJl + A [2t + (k - k2)2 + (k3 - k)2] 



icus + A [2t + k2 + (k3 - k)2] -ZCJ2 + A [2t + k2 + (k - k2)2] 



, (84) 



where the short lines indicate the vertices for clarity and do not enter the expressions of the vertex functions, k^ = 
ki + k2, and a;3 — lui + uj2. One can check that when all external frequencies become zero, the vertex functions 
indeed recover their static counterparts. For example, the integral in Eq. (82) remains unchanged after a replacement 
k' — > k — k' . Adding these two integrals up results in 



r„(k, 0) = A(. + k^) - Ixgl J dk'^^^j^, 
which is Ar|' (k) . Similarly, one can show that 

ri2(k„0) = Xg, - Ag|/^k ^^^j,2,+ (i,!-k)2r+(k!+k)2 - 



(85) 



(86) 



which is just Ar|*(ki). Moreover, one can also check that An effect of the interaction is the shift of the instability 

the fluctuation-dissipation theorem Eq. (70) does obey point. It is defined by 
to this order by a similar manipulation. 



rii(o,o) = 0, 
rio(o,o) = A/i„ 



(87a) 
(87b) 



11 



which are 



^ 2 
rs = 2-93 



1 



1 



-.93 



(T.+k2)2' 
1 



(88a) 



dk- 



2"" 7 """Ts+k2 

for a uniform external field H and hence h and M. 
Eqs. (87a) and (88a) determine the point at which the 
static susceptibility x diverges (see below). Eqs. (87b) 
and (88b) specifies the field at Tg for 



m{Ts,hs) = {ip{Ts,hs)) = 0. 



(89) 



It can be shown explicitly that Eq. (88) is exactly the 
instability point of the one-loop (p^ free energy at m = 
but the corresponding (j)'^ free energy has a vanishingly 
small barrier in accordance with what observed near the 
spinodal point in non-mean-field systems. Accordingly, 
this instability point has been advocated as a new defi- 
nition of a general spinodal point [19]. As the mean-field 
spinodal given by Eq. (39) is Ts = 0, one can replace r, 
in the integrand in Eq. (88a) by its mean-field value in a 
perturbation expansion. So, 



1 



dk 



1 

k4' 



2.93 



^4- 



(90a) 
(90b) 



The first expression is infrared divergent for c? < 4 and 
the second for d < 2 because they diverge when |k| — > 
in these dimensions. Note that the momentum integral 
is cutoff at A. 

The infrared divergences in Eq. (90) affect the posi- 
tion of the instability point which is a property of a spe- 
cific system. They can in fact be subtracted by mass 
renormalization as we shall do shortly. However, these 
infrared divergences which stem from the small momen- 
tum behavior at the instability point persist and in fact 
plague the whole perturbation expansion below an upper 
critical dimensions. 

Equation (90a) can be employed to write the static 
susceptibility x as 



X 



rii(O,0)/A = r, 



1 1 2 



dk 



2k2 



(T^+k2)2k4 



since the correction incurred is of higher order in 93, 
where Tr = t — Tg. Note that the term proportional 
to Tr in the integral equals zero when = 0, one sees 
then one can again set in the integrand to its mean- 
field value if c? > 6 because the integral is then finite. As 
a result, x is again proportional to t~'^ with the mean- 
field result 7 = 1. However, if d < 6, the integral has 
an infrared divergence at the instability point and the 
mean-field result becomes problematic. In fact, in this 
case, the integral is dominated by small |k|s and so can 
be well approximated by 



1 



r(3 



(47r)'i/2(d-4) 



(92) 



by extending the momentum cutoff to infinity and using 
Eqs. (Bll), (B2), and (B3), where T here is the Euler 
Gamma function. One sees from Eq. (92) that the second 
term in the brackets can become arbitrarily larger than 
1, the first term, for sufficiently small t,. and the per- 
turbation expansion breaks down. This therefore iden- 
tifies the upper critical dimension dc — Q above which 
the mean-field theory of the ip^ model is valid. Accord- 
ingly, the spinodal curve and also the mean-field spinodal 
fixed point [41] which build on the mean-field theory all 
become valid above dc in agreement with previous mean- 
field analyses [41, 42]. 

Moreover, one can even introduce a Ginzburg effective 
temperature tq to write Eq. (92) as 



Tr 



(93) 



similar to the critical phenomena [60, 61]. Accordingly, 
for d > 6, sufficiently small r^, i.e., sufficiently close to 
the instability point can always make the second term 
in the square brackets smaller than the first one and the 
mean-field theory is valid. However, for d < 6, this is 
only true for jT^j ^ tg, which is the Ginzburg criterion 
[i't'l]. This therefore defines an unstable region similar 
to the critical region within which fluctuations dominate 
and the mean-field theory fails. Note that, from Eq. (92), 
this is always true in d = 4. This special dimension would 
have become 2 [()2] if the term proportional to in the 
integral in Eq. (91) had been neglected. 

We shall show in the following (Sec. VII B) that higher 
order corrections also confirm this differentiation of spa- 
tial dimensions. Note that our analysis here is completely 
parallel to that in critical phenomena [45, (i3, 64]. There- 
fore, the natural solution to the infrared divergence prob- 
lem below dc is the RG theory. 



VII. RENORMALIZATION-GROUP THEORY 

We have shown in the previous sections that an FOPT 
can be mapped to the instability point of a (p^ model 
and the long wavelength fluctuations near the point re- 
sult in strong infrared divergences below dc in close anal- 
ogy to the critical phenomena. In this section, we shall 
apply the RG theory in the field-theoretic formulation 
[43-47, 63, 64] to the (p^ theory. The RG technique is an 
effective method for treating both infrared and ultravi- 
olet singularities [4(3]. The usual field-theoretic formula- 
tion often concerns with the latter. We shall discuss their 
relation by applying the usual approach to the present 
theory and formulate the RG theory for the transition 
depicted in Fig. 1(c), an RG for a massless theory, i.e., 
the RG for the pfl theory at the point where Eq. (87a) is 
satisfied similar to the critical massless theory. 
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A. Canonical dimensions 

To start with, we determine the canonical dimensions 
of various quantities [45, 64]. These can be obtained 
by a naive dimensional analysis in which 'H, "Ha, and £ 
arc regarded as pure numbers as they appear in expo- 
nentials and all other quantities have the dimension of 
length, which is taken to be —1, i.e., [|x|] = —1, where 
the square brackets here denote the canonical dimension 
of the quantity in them. So, [|k|] = 1. Eq. (58) then 
leads to [Xt] = —2 and [r] = 2 as d/dt/X and r play a 
similar role to V^. Consequently, [(f] = (d + 2)/2 and 
[if] = {d-2)/2 from Eq. (58) and hence [A J] = {d + 6)/2 
and [A J] = {d + 2)/2 from Eq. (21), where the A factors 
arise from the time integral. As a result, [53] = {Q — d)/2, 
[h] = [XJ] = {d + 2)/2 from Eqs. (17) and (8). Simi- 
larly, one finds [r] = 2, [0] = (d - 2)/2, [$\ = (d + 2)/2, 
[H] = {d+ 2)/2, and [g] = 4 - d for the dynamic 0"* 
model. Note that [gs] = [g] + [M] since [M] = [cj)]. 
These dimensions are identical with their static coun- 
terparts as can be readily seen from Eqs. (1) and (8). 
In addition, from Eq. (59), [Av3(k,a;)] = -(d + 6)/2 and 
[Xf>{k,uj)] = -{d + 2)/2. 

Accordingly, one finds from Eqs. (25), (27), and (30) 
that [Gnn'{{xi})] = n[(j)]+n'[(f>] = {n+n')d/2+n' -n and 
[A"'+"r„-„({a;J)] = n[XJ]+n'[XJ] = {n+n')d/2+n'+3n, 
the latter can also be obtained by attaching back the 
amputated legs of Tn>n{{xi}) with integrations and re- 
lating it to its corresponding response function. We 
have used {xi} to represent the set of all the spatial 
and temporal arguments. Then, from the Fourier trans- 
forms defined in Eq. (63) and a similar definition for the 
vertex function, one obtains [A""*"" Gnn/ ({k^}, {wi})] = 
[Gnn'{{x^})] - {u + n'){d + 2) = -{n' + n)d/2 - n' - 
377, = n[A<(2(k,a;)] + n'[Ai^(k, w)] and [r„/„({ki}, {w^})] = 
[A"+"'r„'„({a;J)]-(n-hn')(d-H2) = -{n' +n)d/2+n-n' . 
The response functions and vertex functions we have pri- 
marily used in Sees. V and VI and shall be used below 
have, however, an extracted delta function factor of di- 
mension —d — 2 arising from the global momentum and 
frequency conservations due to space and time transla- 
tion invariances as seen in Eq. (63). Consequently, their 
dimensions change by this amount and are listed in Ta- 
ble II using the same symbols but with one power of A 
less due to the same reason. 



B. Power counting and divergences 

We now consider the relation of the infrared diver- 
gences to the ultraviolet ones and show that the con- 
sequence of the one-loop analysis in the perturbation ex- 
pansion above survives to all orders. In particular, the 
ip^ theory is non-renormalizable and thus the mean-field 
results are valid for d > dc but super-renormalizable for 
d < dc in which the mean-field theory fails and only the 
cubic interaction is relevant if a massless theory exists. 



Take Eq. (88b) as an example. Assume that a mass 
renormalization has been performed such that all rs in 
the response functions have been replaced with t^s, which 
vanish at the instability point. The integral diverges for 
d < 2 because it behaves as |k|'^~^ when = as can 
be readily seen from Eq. (90b). If we change the variable 
to k/-y/7v, we have 

hs = \9zrt'^-' j rfk^^, (94) 

where the cutoff now becomes K/^Jt^. For d > 2, 
the integral is dominated by large |k|s and diverges as 
{A / ^Jt^Y""^ . One sees, however, that in this case the 
Tr factor just cancels the one outside the integral and 
hg ~ A'*^^, a constant. So, no divergence appears as 
should be. However, for d < 2, the integral is now domi- 
nated by small |k|s and can be integrated using Eq. (Bl) 
when T,. — > 0, leading to 

This result coincides with what would arise from a direct 
integration of Eq. (88b) for small and thus describes 
the infrared divergence at the instability point correctly. 
However, as the new cutoff is now infinity at the instabil- 
ity point, the integral can be directly integrated without 
the approximated argument that leads to Eq. (92). Ac- 
cordingly, we can study directly the ultraviolet behavior 
of large A in both cases. This enables us to employ the 
machinery of the field-theoretic RG theory for ultraviolet 
divergences. From the point of view of statistical physics, 
this RG procedure singles those universal properties out 
of the microscopic details specified by A. 

Before extending the analysis to all orders of the per- 
turbation expansion, we first consider in general the 
asymptotic dependence of r„/„({ki}, {uii}) on A, i.e., its 
large momentum behavior. We do this by only consider- 
ing primitive dependence and neglecting subintegrations 
[45, ()4] as an illustration. At the A^th order, a graph 
of L loops contributing to r„'„ contains L integrations, 
N vertices each emitting 3 lines, [3A^ — {n + n')\/2 in- 
ternal lines among which N — n' are response lines each 
having a k~^ factor and the remaining [N ~ n -\- 'n')l2 
correlation lines each having a k^"* factor for large |k|s. 
Consequently, the graph has a total k factor of 

5n-n = L{d + 2) - 2{N - n') - 4 X i(7V - n + n') 

^ L{d + 2) - m + 2n (96) 

powers and thus is proportional to A''"'". Note that the 
frequency integrations are also relevant to momentum 
factors as can be checked in Eqs. (81) to (84). The L in- 
ternal integrations comes from the total internal lines less 
the number of momentum and frequency conservations at 
each vertex except for one for the overall conservation as 
may be seen in Eq. (63). So, 

L = ^[3N~{n + 7i)]-N+l = ^{N-n-n') + l. (97) 
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TABLE II. Canonical dimensions 

|x| |k| \t T v{x)/M g 33 h/H/\J \J \ip{\L,i^) Ay(k,a;) A'^+^'^Gw ({k,}, {a;.}) \-^V,,,,,{{\<i,},{uj,}) = d,,,. 

-11-2 2 ^ A-d^ ^ ^-"^ "" d - n - 3n + 2 ^""^ "" d + n - n' + 2 



Combining Eqs. (96) and (97) leads to 

5n'n n')d + n - n' + 2 +^{d - &)N 



(98) 



where we have use the dimensions of the two quantities 
given in Table II. The last expression is apparent as the 
dependence on A must be deducted from the total di- 
mension the contribution from the coupling constant 173. 
All these results may also be checked with Eqs. (81) to 
(84). 

Equation (98) has one overall factor depending only on 
d and on n and n' of the vertex function and the other 
depending linearly both on the order of the perturbation 
expansion and on the dimension of the coupling constant. 
In fact, as the first factor does not depend on the cou- 
pling, it can be shown to be the dimension of the vertex 
function in the Gaussian theory. One concludes therefore 
that if the coupling constant is dimensionless, the cutoff- 
dependence of the vertex functions will be independent 
of the order in the perturbation theory. So will the re- 
sultant primitive divergences of all the vertex functions 
when A — >■ cxj. This identifies again a critical dimension 
dc = 6 at which the coupling constant is dimensionless 
for the ip^ model considered. For d > dc, the ultraviolet 
divergences of the graphs of a vertex function increase 
their number with the order of perturbation expansion. 
As a result, they cannot be absorbed into a finite num- 
ber of renormalized parameters which contain the strong 
A or system dependence. The theory in such a case is 
then termed non-renormalizable. On the other hand, for 
d < dc, the degree of the primitive divergences decreases 
as one goes to higher orders in the perturbation the- 
ory and the theory is then termed super-renormalizable, 
while a theory in d = dc is renormalizable [45, 63, 64]. 

Now we relate this asymptotic dependence to the in- 
frared behavior of the graph which contributes Dn'n to 
Tn'n- For simplicity let us set all its external momenta 
and frequencies zero and change again its integration 
variables by 



As a result. 



■ Uj/Tr- 



(99) 



(100) 



The remaining integral in -D,J,/,j has no infrared diver- 
gences at Tr — >■ as can be seen from Eq. (94). However, 
when d > dc, i.e., the interaction is non-renormalizable, 
for sufficiently high orders in the perturbation expan- 
sion, D' , has ultraviolet divergences proportional to 



[A/^Tv]*"'". Yet, similar to our example, the two fac- 
tors of Tr cancel exactly, leaving only finite numbers in 
every order in the perturbation. As a result, the mean- 
field behavior survives. For those vertex functions with 
large n' and n (for d > 2), there may exist some lower 
order diagrams that have Sn'n < from Eq. (98) and 
are thus ultraviolet convergent. As Sn'n increases with 
the order N, the leading infrared divergence comes from 
the lowest one-loop term with only vertices having the 
smallest c?„/„ [45]. 

When d < dc on the other hand, the theory is super- 
renormalizable and only a finite number of diagrams have 
ultraviolet divergences. These divergences can then be 
incorporated into a redefinition of the model parameters 
such that the theory built on these redefined parameters 
converges at large momenta. This means that the large 
scale behavior of the theory does not depend on its mi- 
croscopic details prescribed by A and thus universality 
ensues. 

Note that, for d < dc, except those ultraviolet diver- 
gent diagrams, the others are superficially convergent by 
power counting and the infrared divergences are thus 
given by r^"'"^^ or may be even more singular when 
Tr — ?> 0. Moreover, as Eq. (98) shows, these singularities 
increase without bound with the order of perturbation. 
However, we now demonstrate that the most infrared di- 
vergent terms order by order in the perturbation theory 
are described by the ip^ theory [45, 03, 64]. All higher 
order couplings are irrelevant as far as the infrared be- 
havior is concerned. 

Consider, for example, a theory that contains l types 
of interaction vertices each emitting one (f line and — 1 
(f lines. For such a theory, Eqs. (96) and (97) change to 

5'n,„ = L'{d + 2) - 2 ^ N,n, + 2N + 2n, (101) 



L' = - 
2 



N,n, - {n + n') 



-iV+1, (102) 



and so. 



Sn'n^2^2-n-n')d + n-n' + 2 



+ ^ N, (^n^d - Til - d 

= dn'n - ^TVjfft], 



(103) 



where N = ^^.d [g^] = + d — n^d/2 is the dimen- 

sion of the coupling with fields with = g. Now 
one sees that — [(?t] = n^^d — 2)/2 — d increases with 
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rit for d > 2. The minimum thus gives minimum 
n^{d — 2)/2 — d and hence minimum which is the 

most negative value and thus characterizes the most in- 
frared divergent diagrams in each order of the perturba- 
tion expansion. In the case of the critical phenomena 
of the Ising model, its order parameter has an inversion 
symmetry. The smallest is thus 4 and thus the (j)* 
model reproduces the sum of the most infrared divergent 
contributions order by order in a mean-field expansion 
below its critical dimension dc = 2nJ{'n^ — 2) = 4 at 
which [g^] = [45, 63]. In systems with FOPTs, as their 
symmetries have already been broken, the smallest 
is thus 3 and thus all other terms including terms with 
more spatial derivatives which give rise to k factors in the 
numerator of the momentum integrals are infrared irrele- 
vant for d < Q. At d = 6, mean-field behavior is modified 
by logarithmic corrections as in critical phenomena. 



C. RG scheme and renormalization constants 

We have seen that for d < 6 the infrared singularities at 
the instability points of an FOPT can also be described 
by a (p^ theory. In these dimensions, the ultraviolet di- 
vergences of the theory are super-renormalizable and can 
thus be absorbed in a few parameters by renormaliza- 
tion. There are usually two perturbation RG schemes 
to deal with these divergences. One uses the e = 6 — d 
expansion and the other works directly in a fixed dimen- 
sion. We shall adopt the first one. We shall only study 
the massless theory as mentioned. This is because we 
have been focusing on the FOPT described by Fig. 1(c). 
Although the massless theory has problems associated 
with infrared divergences and cannot yield results such 
as universal amplitude ratios [(iG-70], it can neverthe- 
less give exponents that we concern here. We shall apply 
the technique of dimensional regulation to the ultraviolet 
divergent integrals and utilize the minimal RG method 
to subtract and just subtract the resultant dimensional 
poles [71]. In fact, this technique of analytically contin- 
uation in the number of space dimensions and pushing 
A to infinity has been argued to be a more natural ap- 
proach to the critical phenomena and similarly to the 
'unstable' phenomena here as one can then focus on the 
infrared behavior alone without resorting to the ultravi- 
olet divergences in a theory with a finite cutoff [•'34] . The 
RG method of minimal subtraction has an advantage 
of decoupling dynamics from statics without having to 
choose deliberately renormalization conditions [-VJ]. As 
a result, the static renormalization factors can be chosen 
to be identical to those of the corresponding equilibrium 
model. This can be used as a self-check and can also 
simplify manipulations. 

As we shall use the e expansion, we shall renormalize 
the theory at dc = 6. The first task is to identify those 
vertex functions that are superficially divergent. They 
are given by 5„'„ > at dc, which, at this dimension, is 
just dn'n of the pertinent vertex functions. From Table II, 



one finds these vertex functions are Fig, which is quartic 
divergent, Fn, which is quadratic divergent, and F20 and 
Fi2; both of which arc logarithmic divergent. So, one 
should in principle introduce four renormalization con- 
stants to absorb the divergences. However, they are not 
all independent. F20 is related to Fn by Eq. (70). Con- 
sequently, we need introduce only three renormalization 
constants. Of course, we have to perform subtractions 
or mass renormalizations to deal with the quadratic and 
quartic divergences of Fig and Fn. However, as we shall 
study the massless theory and utilize the technique of di- 
mensional regulation, these subtractions become zero be- 
cause massless integrals such as Eq. (90) are zero in the 
dimensional regulation. This does not mean, of course, 
that the instability point is not shifted. Note that we 
shall not consider composite operators and the free en- 
ergy. 

Although we shall work with the massless theory, we 
shall consider an external field that drives the transition 
as Fig. 1(c) shows. Accordingly, we rewrite our relevant 
action as 



£ = / dxdtip 



dt 



1 2 
2.93^ 



h 



(104) 

where we have neglected and hg as they are zero in 
the dimensional regulation. The rcnormalizability of this 
massless theory for d < 6 means that by choosing the 
renormalization factors Z{uii,e)s as 



ryl/2 



U = ZuUr, 



X = ZxX 



(105) 



the renormalized G„„' n and r„'„fl as functions of the 
renormalized functions and ur satisfying 



Gc 



ryn/2 ryn 



(106) 
(107) 



from their definitions arc finite at every order in an ex- 
pansion in ub. as e — >■ 0, where is an arbitrary momen- 
tum scale, TVd is given by Eq. (B2), and the subscripts 
i? denote renormalized quantities. As the fluctuation- 
dissipation theorem must hold in both the unrenormal- 
izcd and the renormalized forms, one finds from Eqs. (70), 
(105), and (107) 



Z. = Z^PZ-}!^. 



(108) 



This confirms that only three Z& arc needed. 

We now compute these factors to one-loop order using 
the method of minimal subtraction and show that they 
indeed make the vertex functions finite. 

Firstly, we choose Eq. (85) to fixed Z^. Note that all 
T must be set to zero when we refer to the equations in 
Sec. VI in the following sections as we are considering a 
massless theory. Using Eq. (B12) and noting that F(2 — 
d/2) contributes an e pole from Eq. (B4) and hence the 
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other r functions can be evaluated directly at d = 6, one 
finds 



rn(k,0) = Ak2 + lA53'A^d(k2)' ' [1 + 0(e)] . (109) 



As mentioned in Sec. VI, rii(k, 0) = Ar|'(k). So, using 
the definition of u in Eq. (105), we have 

rr(k) = k2|i + lz.Mi + o(e)]|, (110) 

where we have expanded (Ikl///)^"^ = 1 — eln(|k|//i) + 
O(e^) and neglected all the terms that add to order 0{e) 
in Eq. (110) consistently. One sees that the neglected 
ln(|k|//i) can never overpower the overall k^ factored out 
in Eq. (110) in the infrared limit and thus one has a 
rcnormalizcd massless theory in every order in u such 
that rii(0, 0) = 0. One sees also that the residue of the 
pole is independent of the external momentum. In fact, 
this is a general feature by which the residue of the pole 
of highest order in a graph must be independent of the 
external momenta and the momentum dependence of the 
poles of lower order must be canceled for any set of exter- 
nal momenta [G4]. This offers a check to the calculations. 

According to Eq. (107), one can choose in the 
method of minimal subtraction to subtract and just sub- 
tract the pole in Eq. (110) such that 



is finite. This leads to 



Zu, = 1 u 

^ 6e 



(111) 



(112) 



Next, we choose Eq. (82) with k = to fix Z;^. Using 
Eq. (Bll) formally and factoring the e pole out, one finds 



1 f iuj \ ^ 

= zu[l + ^u^ [1 + 0(1)]|, 



(113) 



where the definition of u in Eq. (105) has been used and 
powers of ln(ia;/2A/^^) has also been neglected. One ob- 
serves similar features to the previous case here. 
According to Eq. (107), we then choose 

{Z^Z^f2 = l-^u' 

to subtract the pole in Eq. (113) and obtain 

1 o 



Z{p 



1 



3£ 



(115) 



using Eq. (112). 

Equation (115) can be checked using the usual method. 
It uses 



5rn(k,w) 



diiu 



= 1 



1 



= 1 + -iru 
4e 



[l + 0(e)] (116) 



from Eqs. (B13) and (B21). This leads therefore again 
to Eqs. (114) and (115). One can also use directly 



r2o(k,0) = -2A- A.g3^/2(k) 



1 



= -2A<; 1 + —u^ 

4e 



[1 + 0(e)] I (117) 
from Eqs. (B14) and (B22). So, Eq. (107) results in 



r20i?,(k, 0, /.t) = Z^r2o(k, 0) 



1 



2XbZxZ^ I 1 + — [1 + 0(e)] ^118) 



4e 



with the help of Eq. (105). Therefore, one has 



Zxz^ = 1 - = {z.^z^Y/'^, 



(119) 



which recovers Eq. (108) as it should be. Of course, be- 
cause the latter is true, Eq. (119) gives rise again consis- 
tently to Eq. (115) using Eq. (112). 

Finally, we utilize Eq. (86) to determine Z^. One has 

T%{k,ii) = ri2«(k,o,/i)/Afl - z3/2ri2fl(k,o,M)/A,j 
^zy^ [gz + glhikiM)] 
= zy^Z^unN-^'^'y^'l'' |l + i^^^l + 0(e)]| 

(120) 

from Eqs. (107), (B15), (B23), and (105). So, 



which leads to 



Z3/2z„ = 1 - -V?, 
^ e 



Z^ = \-^u\ 
4e 



(121) 



(122) 



using Eq. (112). 

However, all the Z factors must be series in ur instead 
of u itself. Using Eqs. (105) and (122), one obtains 



4e 



(123) 



(114) by iteration. Therefore, to the one-loop order, all can 
simply be replaced by u|j and we finally have 

2 2 ' / ^ 2 

V 6e ^' ^ 3e ^' " 4e 



from Eqs. (112), (115), and (122). 



(124) 



D. Renormalized massless theory 

We have obtained a finite renormalized massless the- 
ory by incorporating the dimensional poles to the three 
Z factors. We can then in principle determine unstable 
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properties using this theory. For example, one can use 
Eq. (85) to find the instability exponent t] which is just 
defined as r(k, 0) — > |kp~'' for small |k| at the insta- 
bility point as can be seen from Eq. (76). Indeed, using 
Eqs. (107), (105), (124), and (B21), one finds, to one-loop 
order, Eq. (109) becomes 



TnR{k,0,^l) = (Z^Z^)i/2rn(k,0) 



1 

6e 



1 + ] k 



1 + ^4 + 0(4,) 



1- i^ ln 



7 ^ 
12" 
"Ikl 



with 



V=qUr 



(125) 



(126) 



for small fxs in equivalence to small momenta if one re- 
places k by ^k. In Eq. (125), the first line is an example 
of Eq. (107) and the second line agrees with Eq. (Ill) for 
the static vertex function. Accordingly, one sees clearly 
that the result of the fluctuation-dissipation theorem, 
Eq. (108), is in fact to ensure the correct static limit 
at w = 0. The third line in Eq. (125) exhibits no poles 
correctly as they just cancel among themselves. We have 
factored out a non-divergent part here and exponentiated 
the terms in the braces in the last line. The manipulation 
appears somehow of brute force to this order but shall be 
born out by the RG analysis. 

Similarly, one finds from Eqs. (113) and (76) 



2Afl^2 



+ o(4) 



As a result, one has 



l--4ln2 + 0(4) 



z = 2 - 7/ 



12 



(127) 



(128) 



to one-loop order in the small /x limit. 

These results indicate that one should study whether 
there is a fixed independent of fj. in the small limit. 
If u*fi is zero, one recovers the mean-field result. But if 
there is a finite when — > 0, one has finite rj and z 
which are independent of the scale fi and are thus uni- 
versal. 



E. RG equations 

In the massless theory one introduces ^ as a substitute 
for the natural mass scale to define the dimensionless 
coupling constant u. Different values of lead to dif- 
ferent rcnormalized vertex functions. However, they are 



physically equivalent as they are related to each other 
by a finite multiplicative renormalization transformation 
[45, 46, 64]. This can be easily seen from Eq. (107). The 
rcnormalized vertex functions at different fxs are multi- 
plicatively related to the unrenormalized functions of the 
same given bare theory. So, they are also multiplicatively 
related by transformation factors that must be finite be- 
cause the renormalized functions they relate are finite. 
All such transformations form an RG. The infinitesimal 
transformations in the case of the massless theory sat- 
isfy a differential RG equation. It turns out that this 
equation leads to useful results such as scaling and uni- 
versality which lie beyond the perturbation theory. 

Now we set up the RG equation. We first consider the 
case in which h = 0, i.e., exactly at the instability point. 
Since the bare functions are independent of /i, one derives 
the RG equation for the rcnormalized Tn'nR as 



d ^ d 



d 
duR 



1 



1 



-n'y^ - -n'j^ Tn'nR = 
(129) 

by differentiating Eq. (107), where all quantities are 
rcnormalized and so finite and the Wilson functions are 
defined as derivatives at constant bare parameters (tp, ip, 
A, 33) 



9/i 



P{ur) = fJ. 



duR 



, , 9 In Zu, 
1<p[ur) = /i- 



. . 9 In Zu, 
1(p[ur) = n- 



Taking logarithms and then partial derivatives of the 
equations associated with A and u in Eq. (105) at con- 
stant bare parameters and taking into account Eq. (108), 
one finds 

dlnZx 1 1 



dlnZu 



-^SUR 



"fuUR. (132) 



Next, we consider the case in which a uniform external 
field Hr ^ hs+rM + from Eq. (6), i.e., Hr is not 

exactly at /is. In this case, m is not zero. So, Eq. (87b) 
becomes 



A_Rft,fl.(a;) = rio_R(0, 0; TOfl(w)). 



(133) 



We have included the case in which itir depends on time. 
Expanding the right hand side at the instability point at 
which m{Ts, hg) = and Eq. (87b) holds, one gets 

00 ^ 

\RhR{uj, Xr, niR, ur, /i) = — ri„_R(0, 0; Q)mJ},. 

n=l 

(134) 

Therefore, using Eq. (129) and noting that m = Z^ ttir 
similar to (p, one obtains 



d ,9^91 91' 

^^a;. + ^^^^ 9A^ + ^d^- 2^-"^« 9;^ - 



hR 
(135) 
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which can also be directly verified by substituting 
Eq. (133) in it. 

Another method to derive Eq. (135) is to note that 
Eq. (133) can also be written in the bare form with a 



F. Solutions to RG equations, fixed points, and 
instability exponents 

Equations (129) and (135) can be solved by the method 
of characteristics [45, G4]. We consider first the first cqua- 



bare field h. Its transformation can then be found using tion. Along a characteristic parameterized by k &nA de- 

termined by flow equations and their respective initial 
conditions 



Eqs. (105) and (107) to be 



(136) 



K— = n 

ClK 

dK 

dUR{K) 



which assures Xhip — XRhnipn and so XrHr is only a shift 
of the source J. Equation (136) is identical with its static 
form. Taking this for granted and using Eq. (107) for Fio 
again result in Eq. (108). Differentiating Eq. (136) with 
respect to /i then leads directly to Eq. (135). 

I 



dn 

Tn'n satisfies 

dVn'nR 

dK 

whose solution is 



= (137a) 
7AAfl(K), Xr{1) = Xr, (137b) 
/3[ur{k)], ur{1)^ur, (137c) 



^"7tp + Ti'^'lv ) r„'„_R, (138) 



1 dx 

r„'„i?X{ki}, {i^i}; Xr{h),ur{k), fin) = r„'„_R({ki}, {w^}; A_R, u_R, fj,) exp <{ - / — [n'y^ix) + n'-/^{x)] } , (139) 



where we have used the solution of Eq. (137a), /x(k) = fiK. 

On the other hand, from the nai've dimensional analysis in Sec. VII A, one has a generalized homogenous relation 



r„'„fl({kj, {uJi}; Xr, ur, /i) = k''"'" AHr„/„i^({ki/K}, {tOi/ Xrk^}; ur, k). 



(140) 



Note the Xr factor from Table II for dimensional reasons. Applying this equation to the left hand side of Eq. (139), 
one finds 

( 1 dx ] 

r„'„ii({AckJ,{wi}; A_R,Mij,/x) = Afl(K)r„/„ii({ki}, {wi/Afl(K)K^};Mfl(K),/z)K''"'" exp j — [nj^ + "'7^] ^ : 

(141) 

where we have relabeled k with Kk. The n dependence here is totally free and can be suppressed. In fact, one can 
simply set k = fj, and ^ itself, the value at k = 1, just 1 at the beginning. 

Equation (141) relates vertex functions at different wavenumbers. In particular, it quantifies how the vertex 
functions change when the wavenumbers are reduced. No simple scaling occurs, however, as the momenta are rescaled. 
The reverse happens when the coupling constant reaches a value such that any further rescaling does not affect it. 
This value is found at the fixed points satisfying 



P{n*R) - 0. 



(142) 



At such a fixed point, ur does not change with k from Eq. (137c). As they depend on k through ur, 7^, 7^, and 7a 
also assume their fixed point values, which will all be marked by stars. As a result, Eq. (137b) leads to 



XrIk) = Xrk''^ 



and, using Eq. (131) for 7^, Eq. (141) becomes 



T,,,MWA^^h^B.u},)^XR^,,.,R{{k,},{uJXR}K-^-^'-;u^ 
which exhibits exact scaling. In particular, 

riii^(Kk, uj; Xr, u^) = Aflrii_R(k, (w/ Xr)k~^~''^ ; w^)k^"'^^ 

I 



(143) 
(144) 
(145) 



using Table II for c?„ 
finds that 



Comparing with Eq. (76), one 



7a 



V = 7^, 
7? = 7? + 2z - 4 



(146a) 

(146b) 
(146e) 



using Eq. (131). Therefore, Eq. (145) becomes 

rii_R(Kk, w; Xr,u*i^) = A_Rriiij(k, (w/Afl)^"^; u^)/*^"**. 

(147) 
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For w = 0, k = 1, and k= |k|, Eq. (147) reads 

TuRiK 0; Xr, u*r) ^ XRVuRih 0; u*R)\k\^-\ (148) 

which is in fact Eq. (125). For k = and k = (w/Ai^)l/^ 
it becomes 

rnfl(0, uj; Xr, uJ,) = XRYnaiO, 1; uJ;)(c^/Ai^)(2-'')/-, 

(149) 

which can be shown to be just Eq. (127). Accordingly, 
one has 

TnR{KO;XR,u*R)^XR (l-^Jg^) |k|'+*, (150) 



Expanding P{ur) in the vicinity of a simple fixed point 



rii_R(0, Lo\ Xr, u*r) = iXr 1 



12 y VA 



(151) 



to one loop using Eqs. (165) and (166) below. 

Similarly, defining one more flow equation associated 
with niR as 

dmR{K) 1 

K ^--jl^TTLR, TO_R(l)=m_R, (152) 

one arrives at the counterpart of Eq. (141) as 

hR{t,XR,mR,UR) = K^'^+^^/^exp |-i 

x/ifl(Afl(K)tK2, m^j(KK-('^-2)/2^ mk('«)) , (153) 

where we have switched to the time domain and used 
Table II. At the fixed point. 



to_r(k) = niRK 2'' 
from Eq. (152) and hence 



(154) 



hR{t, Xr, niR, u*r) = K'^^^''hR{XRtK'' ,mRK '^^^ ,u*r) 

(155) 

with 



P5/u=\y + 2-i^), 



(156) 
(157) 
(158) 

using Eqs. (143) and (146). Choosing k ~ '^i'r^ in 
Eq. (155) leads to a scaling form for the equation of state 
at the instability point 

hR{t,XR,mR) ^ mRf'i{XRtm''^^'^) (159) 

where /a is a scaling function. 

We have solved the RG equations and studied the be- 
havior of a system whose renormalized coupling constant 
lies exactly at the fixed point. We now investigate how 
and when the renormalized coupling constant fiows into 
the fixed point. 



I3{ur) = (3'{uR - u*r), 
one finds from Eq. (137c) 

ur{k) -Ur^ (ur - u^)k'' , 



(160) 



(161) 



where /3' is the derivative of /3 at One sees therefore 
that if /?' > 0, ur{k) will flow to u*r for k — > indepen- 
dent on the initial renormalized coupling constant. On 
the other hand, if /?' < 0, ur{k) will flow to u*j^ for k — > 00 
for arbitrary initial renormalized coupling constant near 
the fixed point. One calls the fixed point in the first case 
infrared stable and that in the second case ultraviolet 
stable. The reason is that one sees from Eq. (144) that, 
as K — )■ 0, one probes smaller and smaller momenta and 
thus the large scale behavior. 

We have computed all necessary Zs in Eq. (124) to 
one-loop order. Accordingly, Eqs. (130) and (132) result 
in 



/3("fl) - 
7v5 = /3(ufl) 





3 3 


-^euR^ 


l< 




1 2 


duR 


6 ^ 




1 2 


duR 


3 ^ 



(162) 
(163) 
(164) 



Eq. (162) has three fixed points. The Gaussian fixed 
point ujj = is infrared stable for e < or d > 6 and 
ultraviolet stable for e > 0. The other two purely imagi- 
nary conjugate fixed points are 



2 

-3" 



(165) 



Although they are imaginary in values, they are infrared 
stable and thus control the large scale behavior for e > 
or c? < 6. Moreover, at these imaginary fixed points, 

V=-\e. z = 2-^e, ^ = 2 + (166) 

from Eqs. (163), (164), (146), and (158), all are real. In 
addition, the first two exponents agree with the results 
found in Eqs. (126) and (128) provided the fixed point 
value in Eq. (165) is substituted. For e = or d = 6 
corresponding to ujj = 0, all the exponents recover their 
Gaussian fixed point values listed in Table I. This fact 
shows again that the fixed points, though imaginary, de- 
scribe correctly the large-scale fluctuations at the insta- 
bility point at which the real mean-fleld transition takes 
place. Moreover, using a flnite-time scaling with Monte 
Carlo RG method, we have found that the imaginary 
fixed point can affect the RG flows in the temperature- 
driven FOPTs of the Potts models in d = 2 [72]. From 
these exponents, others can also be obtained. We shall 
return to the instability exponents in Sec. VIII below. 
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G. Finite-time scaling 

Although the scahng and universahty behavior found 
in the last section is similar to the critical phenomena, 
it may not be easily observable because of the instability 
and the dynamic nature of the transition. An accessi- 
ble method is to employ finite-time scaling which has 
been proved to be effective in the critical phenomena 
[56, 57]. In this method, one varies the external field 
linearly through the transition point to probe the scaling 
behavior. We derive the finite-time scaling form for the 
FOPTs in this section. 

Let the rate R of the linear driving transform as 



one finds 



R — ZrRR, 



(167) 



(168) 



by assuming the driving form, Eq. (48) with the A in- 
cluded, is identical in both the bare and the renormalized 
forms such that 



hR = \RRRt = Zy^/i ^ zy^XRt 



(169) 



using Eq. (136) and Eqs. (105) and (108). Note that 
although such a driving breaks the time translational 
symmetry, it has been shown that this only introduces a 
possible initial slip near criticality arising from nonequi- 
librium initial conditions that induce a new singularity 
[73]. In finite-time scaling, we always start the driving 
far away from the critical/instability point as mentioned 
in Sec. IV C. Consequently, the initial slip should show 
no effects except that one starts driving near the point to 
study it purposely. Without considering the initial slip, 
the driving brings on no new singularity [5.s]. Accord- 
ingly, the three Z factors introduced in Sec. VII C suffice 
for curing all the intrinsic singularities. This is why Zr 
is related to the previous ones as seen in Eq. (168). As 
Xr and Rr are related, we choose Rr as a variable and 
write the RG equation for Hr as 



d_ 




IrRR 



d 



■/3 



d 



ORr duR 2 



1 d 

lipTllR 



dniR 



flR 

(170) 



with 



7r = M- 



dhiRR 



dhiZr 



(171) 



from Eqs. (167), (168), and (130). The solution at the 
fixed point is then 



hR{RR,mR,u*j,) - K^^'^hRiRRH-'-", 
with 



(172) 



TH = ^ - 7; = ^ + ^(rf + 2 - 77) = 2 + p5/v (173) 
by noting that 

dRR{K.) 
dn 



IrRRin), Rr{1) = Rr, 



(174) 



and [R] = [h] - [At] = {d + 6)/2 using Table II. Uses have 
also been made of Eqs. (146) and (156) and the solution 
of Eq. (174) at the fixed point. Equation (173) is just 
Eq. (50) and Eq. (172) is just another form of Eq. (49). 
In fact, one can derive directly an RG equation for tur 
with Hr as its variable [35]. As a result, one finds 



d + 2-1] 
d + 2-ri + 2z' 

d-2 + rj 
d+2-r] + 2z 

whose one-loop e expansions are 

2 1 _ 1 7 

^ 3 



hh = £, 

3 54 ' 



(175) 
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respectively, using Eq. (166), which recover their mean- 
field values for e = 0. These hysteresis exponents were 
found to be comparable with direct numerical solutions 
of Eqs. (1) and (9) [35]. This again supports the relevance 
of the (fi^ theory to the FOPTs. We shall return to this 
comparison in Sec. VIII below. 



H. Leading corrections to scaling 

When the initial ur does not lie at mJj but is near it, 
there are leading corrections to the exact scaling behavior 
as exhibited in Eqs. (144) and (155). 

Consider the external field for an instance and assume 
that /3 is given by Eq. (160). Let 



IiiZ^{ur) 



dx 

W) 



[7^(x)-ry], (177) 



^r{ur) = Aijcxp ■ 



rriRiuR) 
dx 



itirZ-'/^ur), (178) 



[7a(x)-z + 2] , (179) 



(ur - u*jf ) exp I J dx 



(180) 



which are finite renormalizations that eliminate trivial 
deviations from the fixed point theory and correspond 
simply to a change of normalization of the different scal- 
ing variables [45]. Then, we assume 
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Afl, niR, ur, = ^^^{ui^)hR{t, Afl(ufl), muiuu), u^, ^i)U{t, A_r(w/?), mR{uji),u, fi) (181) 
with the boundary condition U{t, Xii{ufi),mfi,{uR),0, = 1 and substitute it into Eq. (135), resulting in 



9 



■ a , a 1 d 

Ofi oXr 2 oniR 2 



dniR 
1 



U{t, XR{uR),mR{uR),u, fi) = 0, 
/ii?(i, >^R{uR),mR{uR), u*R, = 0. 



(182) 
(183) 



These equations can be solved again by the method of characteristics after expanding U in powers of u with the result 



hR{t,XR,mR,UR) = Z-^'''{uR)KP^'''hR{XRtK\mRK~l^'\u*R) 



1 + ^ u^n'P' U{XRtK',mRK-PI'') 



(184) 



where is a scaling function and uses have been made of Eqs. (156) and (157). A finite-time scaling form of Eq. (184) 



IS 



hR{R, XR,mR,UR) = Z~^/\uR)RP'l''-"hR{\,mRR-P'--",u*i^) 



(185) 



with another scaling function f^. 

Note that the exponent of the corrections, /3', depends 
only on the fixed point and is thus universal. To one-loop 
order, 

(3' = e (186) 

from Eqs. (162) and (165). It is real. However, the 
derivative of 7^ with respect to its argument at u^, 

y^ = ^u*R = ±i^^s (187) 

from Eqs. (163) and (165), is imaginary. Expanding /3 
and "f^ to first order in ur — u)^, one finds from Eq. (177) 

i r2~ 

(188) 

which may be complex depending on ur. 



VIII. EXISTENT v?^ THEORIES, YANG-LEE 
EDGE SINGULARITY, AND RESUMMED 
INSTABILITY EXPONENTS 

We have studied in detailed the field-theoretical RG 
theory of the (p^ theory of the FOPTs in the last section. 
In fact, Hamiltonians that possess such a cubic interac- 
tion have been utilized to model many phase transitions. 
Examples include the isotropic to nematic phase tran- 
sition in liquid crystals [74, 75], systems described by 
the Potts model [54, 75-77], in particular the percola- 
tion problem [54, 55, 75] which is the single state Potts 
model [78] , the Edwards- Anderson model of spin glasses 
[79, 80], a lot of nonequilibrium systems [s j ] such as the 



dynamic isotropic percolation and the directed percola- 
tion [82], the Reggeon field theory for high-energy scat- 
tering amplitudes [S3, 84] which falls in the same uni- 
versality class as the directed percolation [85], the Yang- 
Lee edge singularity [3(i] and its related problems such 
as isotropic branched polymers in a good solvent and 
undirected lattice animals, Anderson localization, and 
directed branched polymers and directed lattice animals, 
which relate to the edge singularity in d + 2 dimensions 
[8(i], d + 2 dimensions [87], and d+1 dimensions [Sn], 
respectively, as well as quantum field theory models in 
particle physics [89] . For the case of (p^ models that are 
directly related to the one studied here, they have been 
classified into two categories depending on whether 'un- 
physical' limits such as a particular limit of state num- 
bers and a purely imaginary coupling constant have to 
be taken or not [90]. Field theories in the category with 
real Hamiltonians and without a particular state num- 
ber limit have an unstable ground state and instanton 
solutions such that a perturbation expansion is not an 
adequate approach [91]. Moreover, the 0^ interaction be- 
comes relevant as the space dimensionality lowers from 
six [91]. These are not true for theories in the second cat- 
egory. Thus, the perturbation series has oscillatory terms 
and hence is Borel summable [92, 93] and the quartic in- 
teraction has been found to be irrelevant [92, 94, 95]. The 
Yang-Lee edge singularity belongs to this class and thus 
the e expansion in the RG analysis provides good results 
[!)()]. In the following, we shall show that the dynamics of 
FOPTs near their instability points falls in the same uni- 
versality class to the Yang-Lee edge singularity and shall 
apply its existent results to extend ours for the FOPTs. 
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A. Yang— Lee edge singularity 

According to Yang and Lee [9()], to study the prob- 
lem of phase transitions, it is necessary to study the 
distribution in the complex fugacity plane of the roots 
of the grand partition function. Under a class of gen- 
eral conditions, these roots lie on a circle in the plane 
for the Ising model and its equivalent lattice gas model 
[97]. For the Ising model, the fugacity is proportional to 
exp{—2H) (note that our definition of H has absorbed 
in it the thermal factor) and so the Lee- Yang circle the- 
orem places the zeros at imaginary magnetic fields. For 
T < Tc, there are zeroes at iJ = in the thermodynamic 
limit and the magnetization as a function of H exhibits 
a jump at H = 0; while for T > Tc, there is a gap of 
width 2iHQ(T) within which no zeros exist at all. It was 
found that the distribution of the zeros was singular at 
the edge of the gap [9S] . This Yang-Lee edge singularity 
is described by [3()] 



with 



M - Ma ^ {H - iHoY 



_1 _ d-2 + ri 
^ d + 2-r]' 



(189) 



(190) 



where Mq is the magnetization at the imaginary field 
iHo. 

Consider a continuous spin Ising model in an imaginary 
external field above its critical temperature, viz., Eq. (1) 
with an imaginary H and r > [36, 88]. As there is 
no spontaneous symmetry breaking for r > but only 
the imaginary external field, the shift in Eq. (4) can be 
only to an imaginary M induced by H. One sees then 
from Eqs. (6) and (7) that r is real but both h and 
are imaginary. Therefore, the leading infrared behavior 
of the Yang-Lee edge singularity is governed by Eq. (8) 
with purely imaginary and h [3(i]. A redefinition of 
if to iip [AC)] which is a dummy variable that will be in- 
tegrated out then turns the Hamiltonian back to exactly 
the form of Eq. (8). Therefore, this singularity for T > Tc 
is indeed described by the same model as the FOPTs be- 
low Tc- As a consequence, for d > dc ~ 6, the classi- 
cal mean-field theory results in cr = 1/2 from Eq. (190) 
and Table I. For d < 6, the e = 6 — d expansions of 
the RG functions for the Yang-Lee edge singularity have 
been computed up to three loops from those of the Potts 
model. Their resummed exponents agree impressively 
with those of the high-temperature series analysis [99] 
even down to the exact result of a — —1/2 for the one 
dimensional Ising model where e = 5 [90]. Also the dy- 
namic critical exponent for the Yang-Lee edge singularity 
has been computed up to two-loop order [100]. 

From the fixed point given by Eq. (165), one sees that 
fixed point value of un itself is thus imaginary. This 
again shows that the infrared behavior of the FOPT at its 
instability fixed point is just described by the same field 
theory as that of the Yang-Lee edge singularity albeit in 



opposite temperature ranges. We can then check that all 
our RG functions agree with previous results [54, 90, 100]. 
Moreover, we now employ these results to estimate the 
instability exponents in the following section. 



B. Resummed instability exponents 

The three- and two-loop results for rj and z of the 
Yang-Lee edge singularity are [90, 100] 

1 43 o /16C(3) 8375 \ . ^, a, , , 
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respectively, where C is the Riemann function. The for- 
mer series has been resummed using a [2/1] Padc ap- 
proximant, a [2/1] Pade-Borel method, and a conformal 
mapping technique to provide rather good estimates for 
even d = 1 or e = 5 [90]. The latter has also been re- 
summed by a [2/ 1] Pade approximant with the aid of the 
result of z in c? = given by 

z(d = 0) = 2 - i[2 + ?7(d = 0)], (193) 

in which r]{d = 0) is obtained by a similar approximation 
to two- loop order [100]. Such a [2/1] Pade approximant 
to 77 has been performed using the exact result of 77 = — 1 
in d = 1 for the Yang-Lee edge singularity to two-loop 
order [99]. 

As we now have a three-loop order of ry, we can form a 
[3/2] Pade approximant to 77 using the same exact result 
in d = 1. The result is 



e 1 - 0.908e - 2.2236^ 
'91 - 1.144£- 1.066^2 ■ 



(194) 



We list in Table III the values for various rfs together 
with the averages of those from various other resumma- 
tion methods [90] . One sees that the former is a bit larger 
than the latter for large ds. This is expected as the value 
in d = 1 has been fixed. Accordingly, as our final esti- 
mates, the results from Ref. [90] are kept for d > 3 and 
are averaged with the present one in d = 2, giving rise 
to a closer value to ?7(d = 2) — —0.78(2) from the high- 
temperature series expansion [99]. From these values, we 
can then compute 5 using Eq. (158). In addition, we have 
also collected in Table III the values of v from 

(195) 



d- 2 + 77 

from Eq. (157), since it has been shown that /3 = 1 for 
the scalar Lp'^ model [90, 101]. For the same reason, from 



7 = /3(<5-l), 



(196) 



7 is just given by 5~\. However, the values of v and 7 and 
even 5 look quite strange for lower dimensions. In fact, as 
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TABLE III. Instability exponents. 
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" Mean- field results. 

^ Equation (194) using the [3/2] Padc approximant. 
Exact result [36]. 

Reference [90] (quoted errors reflect the spread in different resummations). 

Final estimates (quoted errors reflect tlie spread in different resummations). 
' Quoted errors reflect the corresponding spreads in "q and/or z in the final estimates, 
s Exact result [86, 100]. 

^ Equation (199) using the [2/1] Fade approximant. 
' Equation (200) using the [1/1] Fade approximant. 
•i Numerical results [35]. 



we are studying the field-driven case, we would consider 
the field instead of the temperature deviations away from 
the instability point. The corresponding relations are 
then 



-7/(9-5 



with 



i///?<5 



d ± 2 - ?7 



j//3S = l-a 



(197) 



(198) 



from Eqs. (156) and (190) and (196). These two exponent 
ratios, which correspond simply to v and 7, respectively, 
in the context of branched polymers [80] and directed 
animals [88], appear normal as seen in Table IIL 

Equation (194) gives a slightly different 77 and hence z 
in d = as compared to z = 1.614 from the two-loop re- 
sult [100]. As a consequence, the [2/1] Pade approximant 
to Eq. (192) then becomes 



e 1 + 1.790e 
18 1 ± 1.515e 



(199) 



with the results for various ds given in Table III. These 
values are also only slightly different from those that ob- 
tained from the two-loop r]{d = 0) [lOO]. 

As mentioned above, the [3/2] Pade approximant to 
T] yields the smallest estimates for d = 5 to 3 and the 
second smallest for d = 2 as can be seen from Table IIL 
We also notice that the direct [2/1] Pade approximant to 
r] [9U] invariably produces the largest estimates given in 



Table III. This suggests us to form a [1/1] Pade approx- 
imant to z, which is 



z = 2- 



1 



18 1 - 0.0153£' 



(200) 



whose results are indeed all larger than those from 
Eq. (199), although the true values may not necessarily 
lie in between them. Nevertheless, we take their average 
as our final estimates for z for d = 5 to 2 to account for 
possible bias due to the two-point Pade approximant as 
seen in the estimates of 77. 

Knowing 77 and z, we can then compute the hysteresis 
exponents from Eqs. (173) and (175) with the results 
given in Table III, where we have also shown their spreads 
due to ry and/or z as has been done for S and v. Note 
that we have computed all other exponents from rj and z 
rather than resummed their respective e series since only 
these two exponents are independent for the ip^ theory. 

In Table III, we have also included the numerical re- 
sults from direct numerical solutions of Eqs. (1) and (9) 
[■-V)]. One sees that the agreement between theoretical 
and numerical results is remarkable. In fact, as pointed 
out above, even the one-loop results agree fortunately 
with these results [.H-^)], though high-order ones without 
resummations do not, similar to the case of critical phe- 
nomena. Of course, higher-order theoretical and further 
numerical results are desirable. Nevertheless, this agree- 
ment confirms again the relevance of the (p^ theory to 
the FOPTs. Note that although we only compare the 
hysteresis exponents in Table III, the static exponents 
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along with the dynamic exponent z can also be estimated 
[72]. As they comprise the hysteresis exponents and arc 
derived directly from the RG theory, they arc more fun- 
damental. Note, however, that they have nothing to do 
with the transition at the usual equilibrium transition 
point. 



IX. DISCUSSIONS 

We have studied in detail the Lp^ theory for the FOPTs. 
The compelling evidences for scaling in driven FOPTs as 
mentioned in Sec. I and demonstrated so far in this paper 
provides strong evidences for the relevance of the theory 
and its infrared stable fixed points to the scaling. We 
have also shown in Sec. VII F and in particular through 
Eq. (161) that the fixed points are indeed reached inde- 
pendent on the initial coupling when k — )■ 0. However, 
as the fixed points are imaginary, one may wonder how 
such imaginary fixed points could be reached from the 
real physical world, or more specifically, how they could 
be reached for an RG flow starting from real physical con- 
ditions and thus questions the extent of the relevancy. A 
possible other interpretation is then that the scaling be- 
havior is only a crossover affected by the tp^ . For the fixed 
points to affect the fiows, however, the latter should fiow 
sufficiently close to the former. Yet, as the fixed points 
are far away from the real physical parameters at least for 
practical es, it would be hardly possible for them to leave 
their trace in real measurements. On the other hand, at 
least in simulations [35] , scaling is easily found in all spa- 
tial dimensions without any detectable indication of dif- 
ficulty for lower dimensions in which the fixed points are 
farther off. Moreover, as has been pointed out, the hys- 
teresis exponents found agree well with those of the (p^ 
theory. Therefore, such a crossover is unlikely to be true 
if not exclusive at all. We also note that the crossover 
of the pseudo-critical phenomena [53] cannot either be 
a candidate as even in the mean-field level the theory 
is irrelevant as shown in Sec. IV C, although the hys- 
teresis exponents that would be predicted by the theory, 
nn = 3/5 for d > 4 from Table I and uh = 0.5493(12) in 
d — Z [58], are not far away from those listed in Table III. 

In order to see explicitly how the fixed points are 
reached, we now study the flow equations in detail by re- 
garding them as rigorous equations. Consider the flow of 
the coupling, Eq. (137c), with /3(ufl) given by Eq. (162). 
This equation is solved analytically by 



2eu\ 



3iip (k^ 



1) + 2£K^ 



(201) 



for the initial uu = Ufl(K = 1). One sees that u\{k) 
converges to u*^ and correctly for k — ^ when e is 
positive and negative corresponding to d < 6 and d > 6, 
respectively. The denominator of the solution equals zero 
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FIG. 3. (Color online) RG flows from k = 1 to «; = of the 
coupling for a real [d = 5 (stars)] and a complex (all the others 
with d < 6) initial coupling, m|j = 0.1 and ufj = 0.1 + O.Olz, 
respectively. The flow for d = 7 has it|j = 0.1 -|- 2i. For clarity 
of illustration, we have cut off the large values of the flow 
with a real initial value at d = 5 (stars), (b) and (c) show 
the projections of the flows in (a) to the real and imaginary 
plane, respectively, in the absence of the flow for d — 7. The 
legend applies to all panels. 



at a specific kq satisfying 



3ul 



3ul 



■2e 



or, 



-sgn(e)2/3iij 



(202) 



to 0(1), where sgn is the sign function. Because the run- 
ning variable k and the dimension s are real, only real 
ur gives rise to a real kq. For e < 0, kq > 1. As a result, 
it does not affect the flow from k = 1 to k = and hence 
any initial uj^ will converge continuously to the Gaus- 
sian fixed point as expected, as is shown in Fig. 3(a). 
However, for e > 0, kq < 1 and the flow will blow out 
at Kq for real wjj. So, kq is the branch point of ur{k). 
Yet, after approaching infinity, ur{k) acquires an infinite 
imaginary part which then converges to the fixed point 
values as seen in Fig. 3(b) and (c). It seems that the 
divergence spontaneously triggers an imaginary part in 
uji{k) so that the purely imaginary fixed point can still 
be reached as the solution (201) shows. Conversely, one 
may imagine that in order to acquire an imaginary part 
to reach the imaginary fixed point, the initial real flow 
has to enter a divergent "no-man's" region. Nevertheless, 
the flow trajectory jumps abruptly from the real plane to 
the imaginary one as illustrated in Fig. 3(a). However, 
it becomes continuous once ur possesses an imaginary 
part, no matter how small it is. The only difference for 
different values of the imaginary part is the magnitude 
of the convolutions shown in Fig. 3(a). The smaller the 
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former is, the larger the latter. Also, no qualitative dif- 
ference shows between different spatial dimensions below 
dc — 6; the difference is only the value of Kq and again 
the magnitude of the convolution. We have checked that 
the divergence remains when considering the contribu- 
tion from two-loop order, though in that case one has to 
resort to numerical solutions. As other variables depend 
on uii{k), they also exhibit a similar feature. 

One sees therefore that if the initial coupling possesses 
an imaginary value, the fixed points can be continuously 
reached irrespective of its magnitude, while for purely 
real initial values to converge to the imaginary fixed 
points the RG flows have to diverge at a finite scale in 
order to acquire an imaginary part. At present, we can 
speculate the reason for the imaginary nature of the fixed 
points. As noticed in [-V)], it has been shown that the free 
energy of metastable states is complex [,37, 42]. In fact, 
since the (p^ theory is not bounded from below, an an- 
alytical continuation has to be performed in computing 
the free energy (2). In particular, as the free energy di- 
verges for negative if, the integration in this region has 
to be deformed to the imaginary axis. So, one should set 
ip ^ lip in the integrand. This naturally leads to ig^. 
Consequently, an imaginary fixed point appears to be a 
natural and plausible choice for controlling the local un- 
stable transition near the instability point of an FOPT, 
just counter possibly to the intuition that only real values 
are physics. It can also be imagined that the divergence 
at a finite scale of a real initial coupling should also be 
pertinent to this. Moreover, note that the negative p 
is just the direction of the FOPT instead of the third- 
order transition. However, before a solid solution to the 
divergence and imaginary problems, the present theory 
remains a hypothesis, though we have shown clearly its 
relevance to the FOPTs. 



X. SUMMARY 

We have basically repeated the whole theory of crit- 
ical phenomena in the (p^ theory for FOPTs in an at- 
tempt to stress their similarity and differences. Instead of 
lengthy proofs we have tried to illustrate known general 
rules with direct examples. We have studied in detail the 
mean-field theory, the Gaussian theory, the perturbation 
expansion, and the RG theory. Finite-time scaling and 
leading corrections to scaling have also been considered. 
We have also touched on the Yang-Lee edge singularity 
and employed its results to improve our estimates of the 
instability exponents. The main results are as follows. 

In the mean- field approximation, we have shown that, 
for a scalar cf)'^ model below its Tc, the FOPT at the 
spinodal point driven by an external field is governed by 
the transition at the instability point of the correspond- 
ing p^ model, although the FOPT and the third-order 
transition between the two states described by the (p^ 
theory fall into opposite ground states and appear differ- 
ent. Via finite-time scaling, numerical results and ana- 



lytical solutions indeed confirm the relevance of the in- 
stability exponents instead of the critical exponents. The 
Gaussian theory shows clearly that both the correlation 
length and the correlation time diverge at the instabil- 
ity point similar to the case of the critical phenomena. 
In addition, all the mean-field instability exponents have 
been derived conventionally. The perturbation expan- 
sion around the mean-field theory then demonstrates as 
expected that infrared divergences plague and thus ne- 
cessitate an RG theory for d < 6, while the mean-field 
results survive for c? > 6 and only outside the unstable 
region for d < 6. The power counting analysis shows 
that in the unstable region and for c? < 6, the effective 
local field (p^ theory does reproduce the sum of the most 
divergent contributions order by order in the mean-field 
expansion. Detailed computations of the renormalization 
functions to one-loop order and derivations and solutions 
of the RG equations for a general vertex function and the 
magnetic field then show unambiguously that there is a 
pair of complex-conjugate imaginary fixed points that are 
infrared stable similar to the critical fixed point and thus 
control the large-scale universal behavior. Exact scaling 
forms and scaling laws among the instability exponents 
have been derived and the latter have been computed to 
one-loop order. In six dimensions, these instability expo- 
nents recover the mean-field ones, which shows again that 
the fixed points, albeit imaginary, indeed describes cor- 
rectly the large-scale fluctuations at the instability point 
at which the real mean-field transition takes place. We 
have computed two particular explicit forms of the renor- 
malized massless two-point response function to one-loop 
order. The finite-time scaling form and associated expo- 
nents have also been derived, which serve as an accessible 
method to probe the instability point. We have derived 
the leading corrections to scaling and confirmed that they 
are controlled by a universal exponent. We have also 
shown that the infrared behavior of the ip^ theory of the 
FOPTs falls in the same universality class to the Yang- 
Lee edge singularity though the former is for T < Tc but 
the latter for T > Tc- This implies that the e expansions 
can be trusted and the cf>'^ interaction is irrelevant. More- 
over, the two- and three- loop-order exponents of the edge 
singularity have been employed to estimate the instabil- 
ity exponents. The outcomes agree well with previous 
numerical results, confirming again the relevance of the 
p^ theory to FOPTs. It appears that this relevancy is 
unlikely just a crossover behavior. On the one hand, 
although the RG flow with a purely real initial coupling 
diverges at a finite scale, an imaginary part generates be- 
yond that scale and the imaginary fixed point can then 
be reached. On the other hand, if the initial coupling 
acquires a finite imaginary part of whatever magnitude, 
true asymptotic behavior can be surely established. We 
speculate that the imaginary nature of the fixed points 
may be a natural consequence of the instability of the (p^ 
theory. However, further studies are clearly needed. 

We conclude therefore that the instability fixed points 
with their instability exponents of the ip^ theory are 
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clearly relevant to the scaling and universality behavior 
exhibited by FOPTs near their instability points. Fur- 
ther studies are desirable to dispel possible concern with 
the instability point, to clarify the divergence of the RG 
flow and the imaginary nature of the fixed points, to 
search for new classes, as well as to find experimental 
evidences. 



Appendix A: Supersymmetry action 

Following Ref. [45], another method to deal with the 
Jacobian, Eq. (15), is to use the Gaussian integral of 
anticommuting classical Grassmann variables, 



det £■ = J VcDc cxp {cEc) , 
to write the generating functional as 



(Al) 



Z[J,J] ^ / V(j>V4)VcVdexp 



-C 



dx 



[.J(t> + J0 



(A2) 



with the effective action 
£[0, 0, c, c] = ( dx 

where 



dcj)_ 6n 
dt ^ 56 



CnCn' + c„'C„ = 0, for any n, n' 



(A3) 



(A4) 



with c„ representing c or c and use has been made of the 
local nature of H. 

Introducing two new Grassmann variables 9 and 9 as 
additional coordinates to x and defining a superfield ip 
by 

^(x, 0, 0) = (j){x) + -/XOcix) + V\c{x)0 + X904>{x), (A5) 
one can write the action (A3) in a beautiful form 



d0d0dt 



A 



with the definitions 



D 



D 



d_ 

80 



dt 



(A6) 



(A7) 



by noting that the integrations over and select the 
term proportional to 09 due to their anticommuting char- 
acter similar to Eq. (A4). For example. 



Icix) 



d9d9dm[i^] = A / _ 5(x)^ 

J { d(j)(x) 6(j){x) 

The action (A6) is invariant under a supersymmetry 
transformation 



(A8) 



or, in the component form, 

64> ~ a/Ace, Sc = 0, 

which mix commuting and anticommuting fields, where 
e here is an infinitesimal anticommuting number. This 
supersymmetry gives rise to Ward-Takahashi identi- 
ties which result in the fluctuation-dissipation theorem, 
Eq. (34), when combined with causality [45]. It also en- 
sures that the equal-time correlation functions converge 
at large times to the corresponding static ones [45]. 

Renormalization can be directly performed to the ac- 
tion (A6) [45]. Direct power counting confirms that 
the static and supersymmetry dynamic theories have the 
same upper critical dimension. And the supersymmetry 
then leads to the renormalized action of the form 

Cr[iI:r] = j d0d0dtl^^Zx J dxD^PRD^/jR + 'HBli'B] 

(AlO) 

So, only one new renormalization factor besides the static 
ones needs to be introduced. Equation (AlO) also implies 
that the form of the Langevin equation, Eq. (9), is kept 
after renormalization. 



Appendix B: Useful formulas and relevant results 

A useful dimensionally regulated integral is [71] 
dk 



where 



r(|)r(n-|), „ 



(47r)fr(f) 



(Bl) 



(B2) 



is the surface area of a d-dimensional sphere divided by 
{27t)'^ and the Euler Gamma function r(z) satisfies 



r{z + l) = zT{z) 



(B3) 



and has poles at zero and negative integers, as can be 
seen from its expansion [102] 



i-iy 



+ iji{n + 1) + 0{e) 



(B4) 



for zero (note that 0! = 1) and integer n for small e, 
where 



1 dTiz) 
which, for the integer n, becomes 

n ^ 



z ^ \l z + l J 



(B6) 
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with the Euler constant 7 = — = 0.577. Using 
Eqs. (B3) and (B4) at n = 0, one obtains an expansion 



r(l +e) = eT{e) = 1 - + 0{e^). 
From Eq. (B5), one also has 



(B7) 



1 



= - (7 + 21n2)e] + 0(e2), (B8) 

where r(l/2) = and -01 (1/2) = -7 - 2 In 2. 



A useful formula to calculate graphs is to use the Feyn- 
man parameters to write 

1 _ Tja) -A- , 1. Sjl-x) 

(B9) 

where a = J2'^i. and x — J^^i-- Using Eqs. (B9), (Bl), 
and the integral 







one finds 



(BIO) 



k2"i(r + k2)" 2 



r(f-ni)r(n + ni-i) 



r(n) 



/ 



dk- 



k2"i(p-k)2"2 2 
Again, using Eqs. (B9) and (Bl), one obtains 



1^^ ^^,Y/2-n,-n. £11) r (i - '^i) r (f - ^2) r (^i + - f ) 



r(ni)r(n2)r(d - ni - 712) 



^2(P) 



/l(p) = / dk 



dk- 



1 



k2 [k2 + (p - k)2] 
1 



2 = -iV<iA (p 



2Nd/2-3- 



"k2(p-k)2[k2 + (p-k)2 

hiPi^P-i) = [ dk 



1 

k2(pi-k)2(p2+k)2 "2 

r 



= -7VdC/-6r ( ^ ) r ( 3 



r 1^3- 
r I 3- 



where pi = /iki, p2 = /ik2, and 



As a result, 



A 
B 
C 



dx 



dx 



X2 



(1 + x) 



,Ad-3 ■ 



dy 



{x + y)2 3(1 _y)! 



(1 + x) 



d-3 



dx / dy [xki + ykl - (xki - yk2)^] 



(B16) 

(B17) 
1-3 
(B18) 



1 



(Bll) 
(B12) 

(B13) 

(B14) 
(B15) 



A = 1 + 0(6), i?=i + 0(£), C=- + 0(e). (B20) 
So, 



Jl(p) = -[l + 0(£)]iV4p2) 

/2(p) = ^[l + 0(e)]A^4p') 



d/2-3 
d/2-3 



For a small £ = 6 — d, one finds from Eqs. (B3), (B4), 
and (B7) 



/3(Pi,P2) = -[l + 0{e)]Ndfi 

£ 



d-6 



(B21) 
(B22) 
(B23) 



F 3- 



l-^e + 0(£^) 



(B19) 



which has an e pole. Owing to this pole, one can simply 
set d = 6 in the integrands of A, B, and C as only first- 
order poles appear in the e expansion to one-loop order. 
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